The Nature of Starburst Activity in M 82 
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ABSTRACT 

We present new evolutionary synthesis models of M 82 based mainly on observations consisting 
of near-infrared integral field spectroscopy and mid-infrared spectroscopy. The models incorpo- 
rate stellar evolution, spectral synthesis, and photoionization modeling, and are optimized for 
A = 1 — 45 /im observations of starburst galaxies. The data allow us to model the starburst re- 
gions on scales as small as 25 pc. We investigate the initial mass function (IMF) of the stars and 
constrain quantitatively the spatial and temporal evolution of starburst activity in M 82. We find 
a typical decay timescale for individual burst sites of a few million years. The data are consistent 
with the formation of very massive stars (> 50 — 100 M©) and require a flattening of the starburst 
IMF below a few solar masses assuming a Salpeter slope dN /dm oc to -2 - 35 at higher masses. Our 
results are well matched by a scenario in which the global starburst activity in M 82 occurred in 
two successive episodes each lasting a few million years, peaking about 10 7 yr and 5 x 10 6 yr ago. 
The first episode took place throughout the central regions of M 82 and was particularly intense 
at the nucleus while the second episode occurred predominantly in a circumnuclear ring and along 
the stellar bar. We interpret this sequence as resulting from the gravitational interaction between 
M82 and its neighbour M81, and subsequent bar-driven evolution. The short burst duration on 
all spatial scales indicates strong negative feedback effects of starburst activity, both locally and 
globally. Simple energetics considerations suggest the collective mechanical energy released by 
massive stars was able to rapidly inhibit star formation after the onset of each episode. 



Subject headings: galaxies: evolution 
galaxies — stars: formation 

INTRODUCTION 



galaxies: individual (M 82) — galaxies: starburst — infrared: 



In recent years, it has become clear that star- 
burst galaxies are important constituents of the 
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Universe at all accessible redshifts (e.g., Heckman 
1998; Adelberger & Steidel 2000; Giavalisco 2002; 
Shapley et al. 2001, 2003; Chapman et al. 2003). 
However, despite extensive studies in the past two 
decades, a detailed and quantitative understand- 
ing of the starburst phenomenon is still lacking. 
Crucial issues that remain open include the evo- 
lution and feedback effects of starburst activity, 
its triggering and quenching mechanisms, and the 
mass distribution of the stars formed in starbursts. 
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Progress has been hindered by the scarcity of spa- 
tially resolved data. Furthermore, high resolution 
optical and ultraviolet studies are often hampered 
by severe dust obscuration. 

In this context, we have undertaken a study 
of the archetypal starburst galaxy M 82 based on 
near-infrared (near-IR) integral field spectroscopy 
and mid-infrared (mid-IR) spectroscopy. Our ob- 
servations provide detailed information on key fea- 
tures tracing various components of the interstel- 
lar medium and of the stellar population, deep 
into the obscured star-forming regions. The data, 
along with results of nebular analysis and popula- 
tion synthesis, are described by Forster Schreiber 
ct al. (2001, hereafter paper I), laying the obser- 
vational basis essential for the starburst modeling 
presented in this paper. In a closely related study, 
Thornley et al. (2000) modeled the mid-IR line 
emission of 27 starburst galaxies to examine the 
issues of massive star formation and evolution in 
starburst environments. Their results are partic- 
ularly relevant to our modeling of M82. 

By virtue of its proximity and brightness, M 82 
is an ideal target for investigations of the star- 
burst phenomenon. Following the seminal paper 
by Rieke et al. (1980), several authors have applied 
evolutionary synthesis models to M82 (Bcrnlohr 
1992; Rieke et al. 1993; Doane & Mathews 1993; 
Satyapal et al. 1997; Colbert et al. 1999, among 
others). However, a number of issues remain con- 
troversial. For instance, it is still debated whether 
the initial mass function (IMF) is biased against 
low-mass stars (see in particular Rieke et al. 1993 
and Satyapal et al. 1997). It has also been sug- 
gested that the formation of stars with masses 
> 30 M Q may be suppressed (e.g., Puxley et 
al. 1989). Furthermore, in analyzing the relative 
distributions of various gaseous and stellar com- 
ponents, some authors have proposed inside-out 
propagation of starburst activity in M82 (e.g., 
McLeod et al. 1993; Satyapal et al. 1997) while 
others have argued in favour of the opposite sce- 
nario (e.g., Shen & Lo 1995; de Grijs et al. 2000; 
de Grijs, O'Connell, & Gallagher 2001). 

The previous models of M 82 were either opti- 
mized to best reproduce the global properties or 
focussed on selected stellar clusters. Given the 
complexity of M82, a full understanding requires 
spatially detailed information of various diagnos- 
tics tracing different star formation epochs. Our 



new IR data, complemented with results at other 
wavelengths from the literature, allow us to probe 
the star formation history in the central regions 
of M 82 "continuously" in space on scales as small 
as 25 pc and in time up to at least ~ 50 Myr 
ago. With the application of starburst models we 
have developed, we use these data to re-examine 
the issue of the IMF and constrain the spatial and 
temporal evolution of the star formation activity. 
Our model results provide quantitative constraints 
for the triggering mechanisms and feedback effects 
of starburst activity in M82. 

The paper is organized as follows. Section 2 
summarizes the observational constraints. Sec- 
tion 3 describes the starburst models which are ap- 
plied to selected regions in § 4. Section 5 presents 
the spatially detailed modeling. Section 6 dis- 
cusses the results and § 7 summarizes the paper. 

2. OBSERVATIONAL CONSTRAINTS 

Our modeling of M 82 uses the observations pre- 
sented and analyzed in paper I. The data con- 
sist of near-IR H- and if-band integral field spec- 
troscopy at R = A/AA ~ 1000 obtained with the 
Max-Planck-Institut fur extraterrestrische Physik 
(MPE) 3D instrument (Weitzel et al. 1996) and of 
A = 2.4 - 45 fim spectroscopy at R ~ 500 - 2000 
from the Short Wavelength Spectrometer (SWS; 
de Graauw et al. 1996) on board the Infrared Space 
Observatory (ISO; Kessler et al. 1996). The 3D 
field of view is nearly parallel to the galactic plane, 
includes the nucleus, and extends to the west out 
to about 200 pc. The SWS apertures cover re- 
gions including and centered on the 3D field of 
view. The 3D data have a spatial resolution of 
1.5", corresponding to about 25 pc at the distance 
of M82 (3.3 Mpc; Freedman & Madore 1988). We 
complemented these data with results from the lit- 
erature, providing further essential constraints. 

We selected five representative regions for a de- 
tailed analysis. These include the "starburst core" 
of M 82 corresponding to a 500 pc-diameter region 
centered on the nucleus, the 3D field of view cover- 
ing the most intense starburst regions, the central 
35 pc at the nucleus, and two 35 x 35 pc regions po- 
sitioned at the brightest Br7 sources observed with 
3D, which we designated M82:Brl and M82:Br2 in 
paper I (« 10" and 5" southwest of the nucleus, 
respectively). For simplicity, we refer to the latter 
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three regions as the "nucleus," "Bl," and "B2." 

For each region, Table 1 lists the main con- 
straints derived in paper I. These include the in- 
trinsic JT-band and Lyman continuum luminosi- 
ties (Lk and L^ yc ) 3 , the total bolometric lumi- 
nosities from OB and cool evolved stars as well 
as the separate contribution of OB stars (Lboi 
and L°®), the stellar masses (M*), the supernova 
rates (z^sn), the CO 1.62 and 2.29 fim bandhcad 
equivalent widths (Wi.62 and ^2.29)5 and the neon 
fine-structure line ratio [Nc III] 15.6 /im / [Nc II] 
12.8 /im. The values of the L k / '-^Lyo iboi/^Lycj 
L^/L Lyc , M* I Lk, and ^sN/iboi are also given. 
Except for the starburst core, the neon ratio is an 
"equivalent" ratio as discussed below. 

The if -band continuum luminosities Lk and 
the equivalent widths W2.29 are corrected for the 
contribution and dilution by hot dust emission 
(which is negligible for Wx.^)- We assumed that 
70% of the OB stars' bolometric luminosity is de- 
tected at IR wavelengths (A = 5 — 300 /xm) as 
reprocessed emission by dust, and that 30% es- 
capes perpendicular to the galactic disk. The stel- 
lar masses M* for the nucleus and the starburst 
core were computed by subtracting the mass of the 
gas component (mainly in the form of H2) from 
the total dynamical mass. Estimates of supernova 
explosion rates i>sn suffer from rather large uncer- 
tainties (see paper I), so related constraints will be 
used as consistency arguments. Among the nebu- 
lar emission line ratios sensitive to the ionizing OB 
stars available from our data sets, we chose [Ne III] 
15.6/im / [Ne II] 12.8 /xm for this work. For the 
starburst core, the ratio is that measured with the 
SWS, corrected for extinction. For the other re- 
gions, we applied the single-star photoionization 
models of paper I to derive "equivalent" neon ra- 
tios, which are those predicted for the OB stars 
effective temperatures inferred from the observed 
near-IR H to He recombination line ratios. 

We also modeled individual regions correspond- 
ing to rebinncd 1" x 1" pixels covering the 3D field 
of view in its entirety. The set of constraints for 
each pixel, from the 3D data and the 12.4 itm map 

3 We use the following definitions: Lj{ is the luminosity 
in the K bandpass A = 1.9 — 2.5 /im (Wamsteker 1981) 
expressed in units of the total solar luminosity Lq = 
3.85 X 10 26 W while LLyc is the Lyman continuum pho- 
ton emission rate times an average ionizing photon energy 
of 15 eV. 



of Telesco & Gezari (1992), includes the intrinsic 
Lk, L\, yc , L°^, W2.29, and the equivalent neon 
ratio. These constraints are the most useful and 
relevant ones for the spatially-detailed modeling 
based on the results of the selected regions. 

3. STARBURST MODELS 

Our starburst models combine evolutionary 
synthesis of stellar clusters with photoionization 
modeling of the surrounding gas. They are opti- 
mized for applications to A = 1 — 45 /im obser- 
vations of starburst galaxies. They assume that 
a given set of constraints applies to a "homo- 
geneous" stellar population, i.e. either a single 
cluster or an ensemble of clusters with identical 
IMF shape and cutoffs (see Thornley et al. 2000). 

3.1. Prediction of Stellar Properties 

We computed the integrated properties of stel- 
lar clusters using the evolutionary synthesis code 
STARS (Sternberg 1998; Thornley et al. 2000). 
STARS is similar to other codes which have been 
presented in the literature (e.g., Bruzual & Char- 
lot 1993; Cervino & Mas-Hesse 1994; Fioc & 
Rocca-Volmerange 1997; Leitherer et al. 1999). 
STARS performs conventional synthesis, assign- 
ing fixed evolutionary tracks to an appropriate 
set of mass bins. STARS employs the Geneva 
stellar tracks (Schaller et al. 1992); for our anal- 
ysis, we selected the solar-metallicity tracks with 
normal mass-loss rates. Although the tracks for 
enhanced mass-loss rates may provide better fits 
to observations (see Meynet et al. 1994, and refer- 
ences therein), the choice has little impact for this 
work. The differences in the synthctized proper- 
ties of interest here are much smaller than those 
produced, e.g., by varying the IMF parameters 
and would imply changes < 1 Myr in the derived 
ages, comparable to those from the measurement 
uncertainties. We adopted a time-independent 
power-law IMF dN/dm oc m~ a between lower 
and upper mass limits m\ ow and m up . The stars 
are assumed to form at a rate R{%) — Rq e~* b /* sc , 
where i?o is the initial star formation rate (ex- 
pressed in M Q yr~ 1 ), % is the burst age, and t sc 
is the burst decay timescale. 

The integrated spectra of starburst galaxies in 
the 1 — 45 fxm range are dominated at the short- 
wavelength end by the direct light of cool evolved 



3 



stars and at the long-wavelength end by radiation 
of mostly OB stars that has been reprocessed by 
interstellar dust and gas. These populations are 
treated with particular care in STARS. Specifi- 
cally the non local thermodynamical equilibrium 
(non-LTE) model atmospheres of Pauldrach et al. 
(1998), which supersede the earlier non-LTE at- 
mospheres of Scllmaier ct al. (1996), are chosen 
to represent adequately the spectral energy distri- 
bution (SED) of the hottest stars (see also Paul- 
drach et al. 2001, and references therein). Em- 
pirical photometric data are used to better ac- 
count for the properties of the coolest stars. Fi- 
nally, the Geneva tracks for intermediate-mass 
stars have been extended to include the thermally- 
pulsing asymptotic giant branch phase (TP-AGB), 
which can contribute significantly to the inte- 
grated if-band luminosity. STARS also predicts 
the strength of various near-IR stellar absorption 
features. Further details on the above aspects are 
given in appendix A. 

The Pauldrach et al. (1998) models include 
the non-LTE radiative transfer and hydrodynam- 
ical treatment of radiatively-driven, spherically 
expanding steady-state winds and account for 
the combined effects of line blocking and blan- 
keting on the radiative transfer and energy bal- 
ance. They provide excellent matches to observed 
high-resolution far-ultraviolet spectra of hot stars 
(Pauldrach et al. 2001) and nebular photoioniza- 
tion computations which incorporate these model 
atmospheres successfully reproduce the relative 
intensities of infrared emission lines measured in 
Galactic H II regions (Giveon et al. 2002). As 
discussed by Giveon ct al. (2002) and Thornlcy 
et al. (2000), the inclusion of the Pauldrach et al. 
models is critical for the photoionization model- 
ing because their SEDs are significantly harder 
than for hydrostatic plane-parallel LTE models 
such as those of Kurucz (1992). Accounting for 
the TP-AGB phase is most important for the CO 
bandheads because it results notably in larger 
equivalent widths (EWs) at ages > 50 Myr (see, 
e.g., Origlia & Oliva 2000). However, as will be 
seen below, our results do not depend much on 
this feature. 

3.2. Prediction of Nebular Properties 

We modeled the neon ratio produced in nebulae 
excited by the integrated stellar SED computed by 



STARS using the photoionization code CLOUDY 
version C90.05 (Ferland 1996). The nebulae were 
represented as a single gas shell surrounding a cen- 
tral point-like ionizing source. In such "central 
cluster" models, the nebular conditions are speci- 
fied by the gas and dust composition, the hydro- 
gen gas density uh, and the ionization parameter 
U = Qhyc/ (47rr 2 7iHc), where Qhyc is the intrin- 
sic Lyman continuum photon emission rate of the 
source, r is the shell radius, and c is the speed of 
light. 

We adopted the parameter values derived from 
the nebular analysis of paper I: rin = 300 cm~ 3 , 
an effective ionization parameter with log£7 e ff = 
—2.3 dex, and solar photospheric gas-phase abun- 
dances (from Grevesse & Anders 1989; Grevesse 
& Noels 1993, as implemented in CLOUDY v. 
C90.05). As argued in paper I, a more realistic 
representation of the H II regions in M 82 consists 
of a well-mixed distribution of ionizing clusters 
and gas clouds. The average ionization parame- 
ter at the outer boundaries of the gas clouds de- 
rived from the observed properties assuming such 
a "random distribution" model constitutes the ap- 
propriate effective value for modeling in the frame- 
work of the idealized central cluster geometry. 
Interestingly, we found nearly identical logt/ e ff 
throughout the starburst core of M82 on scales 
of a few tens of parsecs to ~ 500 pc. As discussed 
in paper I, this appears to indicate a uniform star 
formation efficiency across the observed regions. 

Following Thornley et al. (2000), we computed 
the neon ratio for bursts with decay timescale 
t sc = 1 Myr by assuming that {/ cff varies pro- 
portionally with QLyc and that the maximum 
value of logt/eff equals —2.3 dcx. For bursts with 
longer decay timescales, we convolved the corre- 
sponding star formation rates R(t\f) with the time- 
dependent neon ratio for the short decay bursts. 
We neglected the effects of dust grains mixed with 
the ionized gas within the nebulae; this has little 
impact on our model results (see § 4.5). 

3.3. Modeling Procedure 

We approached the modeling of M82 by first 
making physical and qualitative arguments to con- 
strain the possible ranges for the parameters, and 
then performed quantitative fits. This led to phys- 
ically plausible and successful models although 
they may not be unique. We first investigated, 
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from the properties of selected regions, the gen- 
eral characteristics of starburst activity: the IMF 
cutoffs and the typical burst decay timescale (§ 4). 
Having determined these parameters and identi- 
fied the most useful age indicators, we then mod- 
eled individual regions across the entire 3D field 
of view to constrain the spatial and temporal evo- 
lution of starburst activity (§ 5). 

Most of our constraints are sensitive to m up , t sc , 
and % , so we investigated these parameters simul- 
taneously. The total stellar mass depends only on 
the low-mass IMF (shape and mi ow ), which we ex- 
amined independently. We used ratios of different 
properties derived from data obtained at similar 
wavelengths along with the CO bandheads EWs, 
minimizing uncertainties from the extinction cor- 
rections and from the distance assumed for M82. 
The exception is ^sN/^boi but we used it mainly 
as consistency check because of the uncertainties 
on ^sn (§ 2). We determined the initial star for- 
mation rates Rq by comparing the predicted lu- 
minosities to the measured values. We kept the 
other parameters constant; the effects of variations 
in the most critical ones will be discussed when 
appropriate. Table 2 summarizes the ranges con- 
sidered or the fixed values adopted for the model 
parameters. 

4. MODELING OF SELECTED REGIONS 
IN M82 

Figure 1 compares the properties of the selected 
regions from Tabic 1 to model predictions. The 
models were computed for m up = 100 M Q and are 
shown for burst decay timescalcs of 1 Myr, 5 Myr, 
20 Myr, and 1 Gyr. The diagnostics considered 
are not affected by the low-mass star population 
so a fiducial m\ ow = 1 M© was here adopted. This 
comparison reveals age differences between the re- 
gions, in particular between the nucleus (older) 
and the regions Bl and B2 (younger). The exact 
ages depend on the upper mass cutoff and burst 
decay timescale, which we constrained first. 

4.1. Upper Mass Cutoff and Burst Timescale 

The neon and iboi/^Lyc ratios are our most 
sensitive diagnostics to the upper mass cutoff. Fig- 
ure 2 shows model computations for these ratios 
with m up = 25, 30, 35, 50, and 100 M Q , and for 
t sc = 1 and 5 Myr. The upper mass cutoff is best 



constrained from the properties of the young mas- 
sive stars only. Therefore, the observed bolomet- 
ric luminosities for selected regions in Figure 2 
exclude the contribution from cool evolved stars 
so that the Lj^/LLyc ratios are compared to the 
models. 

Both diagnostics imply similar dominant OB 
star populations for the various regions consid- 
ered. The small values of the neon ratio suggest a 
relative paucity of ionizing photons with energies 
> 3 Ryd (the ionization potential of Ne + ) while 
the low ^boi/^Lyc ratios imply that lower energy 
Lyman continuum photons are being produced ef- 
ficiently. From the comparison with the models, 
two alternative interpretations are possible: (1) a 
high m up > 50 M Q , a short timescale of at most a 
few million years, and the softening of the ionizing 
radiation field attributable to rapid aging of the 
starburst, or (2) a lower m up down to w 30 Mq, 
with longer timescales possible. Upper mass cut- 
offs below 30 M Q are ruled out, but the data do 
not allow us to constrain unequivocally m up and 
t sc for m up > 30 M©. 

As summarized by Thornlcy et al. (2000), a 
variety of results on local templates of starburst 
regions as well as other large-scale starburst sys- 
tems make high upper mass cutoffs more plausi- 
ble. In M82, the low L^/L Lyc values for Bl and 
B2 can only be reproduced with m up > 50 M©. 
The neon and L®® /-^Lyc ratios measured for all se- 
lected regions lie within the ranges determined for 
the sample of dusty solar-metallicity IR-luminous 
starburst galaxies of Thornley et al. (2000). These 
ranges, indicated in Figure 2, are compatible with 
iriup > 50 M© and differences among the sample 
sources attributable to a range in burst ages. 

In addition, the poor correlation in M82 on 
scales of a few tens of parsecs between the spa- 
tial distributions of the near-IR continuum emis- 
sion and CO bandhead EWs, and that of the ion- 
ized gas (e.g., Br7 and He I 2.06 (im) and IR dust 
emission provides strong evidence for short decay 
timescales in individual burst sites (see the 3D 
images of paper I; also, e.g., Telesco et al. 1991; 
Larkin ct al. 1994; Satyapal et al. 1997). Since 
red supergiants dominate the near-IR continuum 
emission in M82 (§ 4.3), the poor correlation im- 
plies that star formation activity was not main- 
tained locally at high levels longer than the typ- 
ical main-sequence lifetime of the most massive 
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red supergiants progenitors (~ 5 Myr). This situ- 
ation is reminiscent of the nucleus of M 83, where 
Puxley et al. (1997) found that the K-band con- 
tinuum emission and CO bandhead strength are 
spatially anticorrelated with the Eh-7 EW, which 
they interpreted, together with starburst models, 
as indicating an evolutionary sequence and burst 
timescales of 1 — 5 Myr. Massive star formation is 
expected to inhibit further star formation locally 
already after a few million years owing to strong 
stellar winds and supernova explosions disrupting 
the surrounding interstellar medium (ISM). The 
~ 100 pc-size gas hole around the brightest ra- 
dio supernova remnant 41.9 + 58 seems to provide 
a particularly striking example of such effects in 
M82 (e.g., WeiB et al. 1999). 

The most remarkable outcome of our models 
is that for m up > 50 M , the OB star-dominated 
properties of the 3D field of view and the starburst 
core imply short burst decay timescales on large 
spatial scales (~ 100 — 500 pc). This is not due 
to a few regions dominating the integrated prop- 
erties; for instance, Bl and B2 contribute together 
less than 15% to the L\, yc and L®^ measured in 
the 3D field of view. Our results suggest that in 
M82 starburst activity can have strong negative 
feedback effects globally as well. We discuss this 
point further in § 6. 

4.2. Successive Starburst Events 

In the rest of this paper, we will adopt m up = 
100 M Q and a representative t sc = 1 Myr as justi- 
fied above. Examination of Figures 1 and 2 then 
reveals that for each region, there is a systematic 
increase in burst age from that implied by the trac- 
ers of OB stars (neon and L^/L^ yc ratios) to 
that inferred from the tracers of cool evolved stars 
(CO bandhcads). The differences are about 5 Myr. 
The Lk/Lljc ratios, measuring the relative popu- 
lations of these stellar components, correspond to 
intermediate ages. 

Figure 3 illustrates the difficulty of single burst 
models to meet simultaneously the constraints for 
each of the selected regions. The figure compares 
the regions' properties to evolutionary tracks for 
single bursts in diagrams of the neon, Lboi/^Lyc, 
and Lk/L^c ratios versus Wi.iv- At ages corre- 
sponding to the measured values of W2.29, models 
for t sc = 1 Myr and m up = 100 M Q predict neon 
ratios 10 (Bl) to 34 (nucleus) times lower than 



measured, Lboi/ '-^Lyc ratios 6 (nucleus) to 23 (B2) 
times higher, and Lk/L^ yc ratios 4 (nucleus) to 
33 (B2) times higher. The constraints are still 
not satisfied simultaneously when augmenting the 
burst decay timescale to 5 Myr with, in particular, 
predicted neon ratios about an order of magnitude 
higher than measured at ages corresponding to the 
observed values of W2.29- Allowing for lower m up 
down to 50 M Q only worsens the discrepancies. 

The age sequence between the various diag- 
nostics may suggest that the hot massive stars 
and the cool evolved stars belong to populations 
formed in distinct, successive starburst events. In 
this scenario, the "young bursts" account for most 
of the ionizing radiation while the "old bursts" 
dominate the near-IR properties. Rieke et al. 
(1993) first proposed such two-burst models of 
the entire starburst core of M82 for similar rea- 
sons. We thus considered models consisting of 
two sequential bursts, each with t sc = 1 Myr and 
m up = 100 M Q , and separated in time by more 
than one decay timescale. Figure 4 shows similar 
evolutionary tracks as Figure 3 for the best-fit time 
separation and relative burst intensities for each 
region derived in § 4.3. From Figures 3 and 4, it 
is also difficult to disentangle between single and 
two-burst models. However, since the neon and 
^b<d / -^Lyc ratios seem to favour the shortest decay 
timescales, successive short bursts could provide a 
natural explanation for the observed strength of 
the CO bandheads. The two-burst scenario also 
appears plausible in view of the morphology and 
nearly edge-on orientation of M82 (§ 6). 

We remark that given the uncertainties of the 
data and models, it may be possible to find single 
bursts with finely-tuned parameters which repro- 
duce all the observed properties. Exponentially 
decaying bursts with t sc intermediate between 1 
and 5 Myr, for instance, could meet the neon ra- 
tio constraint while violating the Lboi/^Lyc and 
La/L^yc constraints at a level that is more se- 
vere than for the double bursts but probably still 
acceptable. As another example, allowing lower 
m up down to 30 Mq together with longer t sc of 
10 — 20 Myr for single bursts may improve the fit 
to the neon ratios but the mismatch is only re- 
distributed to the other properties. Similar diffi- 
culties in modeling comparable sets of constraints 
have been encoutered for other dusty IR-selected 
starbursts, e.g., by Engelbracht et al. (1998) for 
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NGC253 and Goldader et al. (1997) for a sample 
of (ultra) luminous IR galaxies. With our data and 
models of M82, together with the obvious com- 
plexity of the problem, it is not possible to find a 
unique solution. 

4.3. Burst Age and Strength 

Adopting the two-burst scenario as working hy- 
pothesis, we constrained the age of the young 
bursts from the neon and L^/L^ yc ratios, and of 
the old bursts from the CO bandheads. The initial 
star formation rates Rq (with m\ ow fixed at 1 M ) 
were adjusted so that the total contribution from 
both bursts to £l V c and Lk reproduces the ob- 
served values. For identical IMF parameters and 
burst decay timescales, the values of Rq provide 
a meaningful measure of the relative strengths of 
the bursts. 

For the young bursts, the ages inferred from 
each indicator agree marginally, the neon ratio im- 
plying slightly older ages. This discrepancy would 
not be solved by other choices of m up and t sc . 
The difference is comparable to the burst timescale 
and, in view of the measurements and modeling 
uncertainties, is probably not significant. The 
Lfafi /Lhyc ratio may trace a slightly different pop- 
ulation as it is sensitive to a somewhat lower stel- 
lar mass range than the neon ratio is. The effects 
of a distribution of cluster masses and luminosi- 
ties or of dust grains within the nebulae, neglected 
here, would increase the discrepancy (§ 4.5). We 
therefore assigned equal weights to the diagnostic 
ratios. 

For the old bursts, an accurate age determina- 
tion is only possible for the nucleus. It has very 
deep CO bandheads which the models reproduce 
only at ages of 8 — 15 Myr, consistent with the 
average spectral type K5 I derived from detailed 
spectral synthesis in paper I. The EWs at the nu- 
cleus further support short burst timescales: for 
t sc > 5 Myr, the contribution from older popu- 
lations with shallower bandheads and/or photo- 
spheric and nebular featureless continuum emis- 
sion from OB stars results in predicted EWs lower 
than observed, for any burst age. Figure 5 in- 
dicates that older ages when intermediate-mass 
stars reach the end of the asymptotic giant branch 
(~ 50 Myr — 1 Gyr) are also ruled out from Wi.tz- 

For the other selected regions, we can set firm 



lower limits on the ages but the EWs alone do not 
provide strong constraints on upper limits. For 
comparison with our data and models, we indi- 
cate in Figure 5 the ranges observed by Oliva et al. 
(1995) for old (1 — 10 Gyr) stellar populations char- 
acteristic of the central regions of elliptical galax- 
ies and bulges of spiral galaxies (hereafter "nor- 
mal populations"). The EWs at Bl are somewhat 
lower than for normal populations and may indi- 
cate an early stage when the first red supergiants 
appear. For B2, the 3D field of view, and the 
starburst core, they are slightly higher and imply 
possible upper limits on the age of « 50 Myr and 
on the timescale of ~ 5 Myr. 

The M* /Lk ratio provides an additional con- 
straint. Empirical determinations for normal pop- 
ulations lie in the range 10—30 M Q /L Q (Devereux, 
Becklin, & Scoville 1987; Oliva et al. 1995; Hunt et 
al. 1999). 4 As argued in paper I, the substantially 
lower values ~ 1 Mq/L© measured for M82 be- 
tween radii of « 10 and 250 pc indicate very large 
contributions from red supergiants to the near-IR 
continuum emission, even in regions of low surface 
brightness. Consequently, ages < 50 Myr are the 
most plausible. We adopted the solutions corre- 
sponding to the rising part of the predicted EWs 
with t sc = 1 Myr but older ages of 15 — 30 Myr are 
also possible. Our choice implies lower Rq for the 
old bursts since the predicted Lk peaks around 
10 Myr. The younger age solutions arc unaffected 
by uncertainties related to numerical fluctuations 
due to our conventional synthesis technique and to 
our TP-AGB implementation, which may become 
a concern at ages > 50 Myr (see appendix A). 

4.4. Model Results for Selected Regions 

Tables 3 and 4 report our two-burst model re- 
sults for the selected regions. The derived burst 
strengths outline the particular intensity of star- 
burst activity in the central few tens of parsecs of 
M82 about 10 Myr ago. Our model for the star- 
burst core agrees well with the results of Rieke et 
al. (1993) who found that separate bursts occur- 
ring about 5 and 13 Myr ago, each lasting a few 
million years, could reproduce the observed prop- 

4 Our units for the M* /Lk ratio differ from the conventional 
definition expressed in (M/L^)q. The quoted range has 
been computed in Mq /Lq from the data presented in the 
given references, assuming a K bandpass A = f .9 — 2.5 /im 
and expressing L K in Lq = 3.85 X 10 26 W. 
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crtics. The small differences between our model 
and theirs reflect mainly differences in the adopted 
set of constraints and model assumptions (such as 
their choice of Gaussians to represent the star for- 
mation rate). 

We estimated the formal uncertainties on the 
relative burst ages and strengths i?o by varying 
in turn the observational constraints within the 
ranges corresponding to the measurement errors. 
They are typically < 0.5 Myr and 25% - 40% for 
the ages and strengths, respectively. For the old 
bursts, the minimum age is fairly well constrained 
by the CO bandhcad EWs whereas the maximum 
age from the quoted uncertainties constitutes a 
lower limit since we did not account here for the 
older possible solutions. The absolute burst ages 
and intensities may suffer from larger uncertain- 
ties related, e.g., to the burst decay timescale and 
shape of the star formation rate. The relative ages 
and intensities are less affected because the evolu- 
tion of the properties considered are dominated 
by stellar evolution during the early phases of a 
burst and for decay timescales comparable to the 
main-sequence lifetime of the most massive stars. 

Figure 6 illustrates the results for the most rel- 
evant properties: L K , L Lyc , L hoh W2.29, \V1.e2, 
and the neon ratio. For each region, the curves 
show the variations of the properties as a func- 
tion of the time elapsed since the onset of the old 
burst and are normalized to the observed values. 
They represent the combination of the two bursts, 
weighted by the respective values of Rq for each 
burst. The models for all regions reproduce the 
main constraints within a factor of three or better. 
The masses computed for the nucleus and the star- 
burst core are, however, substantially lower than 
the measured values. The mass is very sensitive 
to the lower cutoff and shape of the IMF, which 
we discuss in § 4.6. 

Our results are roughly consistent with the con- 
straints imposed by 1*3 n and 10 12 i^sN/^boi- A 
good agreement for the regions observed with 3D 
is not expected because of the large uncertainties 
in deriving j^sn from [Fe II] 1.644 fim line fluxes, as 
discussed in paper I. The ages inferred for the old 
bursts satisfy the "timing constraint" implied by 
the extent of the X-ray halo tracing the supernova- 
driven wind of M82, which requires a minimum 
starburst age of about 10 Myr (Heckman, Armus, 
& Miley 1990; Rieke et al. 1993; Lehnert, Heck- 



man, & Weaver 1999). 

4.5. Model Uncertainties 

We investigated the effects of various sources of 
uncertainties on our model results so far, includ- 
ing the assumption of a homogeneous cluster pop- 
ulation, the nebular parameters, the shape of the 
IMF, and the synthetic model atmospheres. Their 
effects have been described for the neon ratio by 
Thornley et al. (2000) whom we refer to for further 
details. Figure 7 shows in addition the effects on 
the iboi/^Lyc an d Lx/L^yc ratios, and on Wi.is- 
None of the above sources of uncertainties affect 
our conclusions concerning the upper mass cutoff 
of the IMF or the ages and relative intensities de- 
rived for the bursts. While possible alternative as- 
sumptions or variations of the nebular parameters 
and IMF shape result in measurable differences 
in the model predictions, the effects remain small 
compared to those of aging or changes in m up . 

In particular, the largest effects are obtained 
for a "heterogeneous" population of stellar clusters 
following the luminosity function (LF) derived by 
Thornley et al. (2000) 5 . This LF was constructed 
assuming it reflects a cluster mass function and im- 
poses a distribution in effective m up of clusters (at 
fixed IMF slope and mi ow ). The large proportion 
of small clusters results in a "down-weighting" ef- 
fect with, for instance, lower neon ratios and shal- 
lower CO bandheads due to the reduction in rel- 
ative population of massive stars. However, the 
effects arc but a factor of two or less. Similar dif- 
ferences for the neon ratio are obtained by varying 
the ionization parameter within the possible range 
-2.5 < logoff < -2 for M82 (paper 1) or adopt- 
ing the softer Kurucz (1992) SEDs for hot stars 
instead of the Pauldrach et al. (1998) SEDs (see 
Thornley et al. 2000). The effects of dust within 
the H II regions, illustrated with the Orion nebular 
composition, are significantly smaller. By absorb- 
ing ionizing photons, dust would also reduce the 
iLyc in the models. This is difficult to quantify, 
but the resulting increase in synthetized iboi/^Lyc 
and Lk /-^Lyc would only imply younger ages and 
keep the inferred burst decay timescales short. 

5 For reference, the LF in terms of cluster H ionizing rate 

Qtyc is dN / d (^Qljc) « {QtycY^ with p = 0.19 for 
Qd yc = 10 45 - lO 49 ' 5 s- 1 and /3 = 1 for Q^ yc = 10 49 5 - 

10 53 B- 1 . 
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Figure 8 shows that the alternative IMFs we 
considered, from Miller & Scalo (1979) and that 
derived by Eisenhauer et al. (1998) for the Galactic 
high-mass star-forming region NGC3603, bracket 
the range of slopes at intermediate and high 
masses determined for young clusters and OB 
associations in the Milky Way and Large Magel- 
lanic Cloud (e.g., Hunter et al. 1997 and references 
therein; Brandl et al. 1996; Massey & Hunter 1998; 
Kroupa 2001). Changes in the IMF shape have the 
smallest effects in our models. As a consequence, 
our set of observational constraints does not allow 
us to investigate the shape of the intermediate- 
and high-mass IMF in M82. 

Potential problems with the stellar evolutionary 
tracks represent an additional source of uncertain- 
ties. A detailed discussion is beyond the scope of 
this paper but we briefly mention two points rele- 
vant to this work. Origlia and co-workers (Origlia 
et al. 1999; Origlia & Oliva 2000) warned that 
model predictions of the CO bandheads can be 
significantly underestimated at ages > 10 Myr 
if the TP-AGB phase is not accounted for (see 
also Chariot & Bruzual 1991; Lancon & Rocca- 
Volmerange 1994) and because of the possible in- 
adequacy of current stellar tracks for the evolu- 
tion of massive stars in the red supergiant phase, 
most severely at sub-solar metallicities. Our star- 
burst models account for the TP-AGB (see ap- 
pendix A) and use solar-metallicity tracks as ap- 
propriate for M82 (paper I), reducing the impact 
of these uncertainties on our results. Furthermore, 
the maximum age from the CO bandheads would 
be most affected, and much less the minimum age 
so that the above uncertainties have little conse- 
quences on the age sequence suggested by our ob- 
servational constraints. This result would be more 
affected if the stellar tracks significantly underpre- 
dict the number of very luminous but relatively 
cool stars (logiboi > 5.75, logT e ff < 4.5), a possi- 
bility raised by Thornley et al. (2000) in their dis- 
cussion of the Galactic Center. With correspond- 
ingly lower predicted neon ratios, the constraints 
might be more clearly accomodated with single 
bursts of somewhat longer durations. However, 
this point remains to be thoroughly investigated. 

4.6. Low-mass End of the IMF 

The low-mass IMF for the starburst popula- 
tion in M82 has long been debated. Various au- 



thors argued that it is deficient in stars with initial 
masses below w 3 M Q compared to the solar neigh- 
bourhood IMF (Rieke et al. 1980, 1993; Bernlohr 
1992; Doane & Mathews 1993). The studies by 
Rieke and coworkers have been the most influen- 
tial. Rieke et al. (1980) argued for a sharp trunca- 
tion of the IMF at 3.5 M©, with no stars formed 
below this cutoff. Rieke et al. (1993) modified this 
conclusion and argued that while the IMF may 
extend down to 0.1 M©, the slope flattens signifi- 
cantly near 3 M©, at a substantially higher mass 
than in the solar neighbourhood (near 0.5 M©; see, 
e.g., Kroupa 2001 and references therein). More 
recently, Satyapal et al. (1997) challenged these 
hypotheses and modeled the integrated if -band 
luminosity of M 82 with a Salpeter IMF down to 
0.1 M© without using up more than w 35% of the 
total dynamical mass in the starburst, a reason- 
able upper limit. 

The origin of the controversy resides primar- 
ily in the notable differences in observational con- 
straints and model assumptions between the var- 
ious studies. Here, we reexamine the issue of the 
low-mass IMF in M 82 in light of our new data and 
models. We consider two regions for this purpose: 
the nucleus (central 35 pc) and the starburst core 
(central 500 pc). 

4-6.1. Reexamination of the Low-mass IMF 

Among our constraints, the stellar mass M* is 
the only one sensitive to the low-mass star popu- 
lation. We computed this critical quantity by sub- 
tracting the mass of the gaseous component M gas 
from the total dynamical mass Md yn (see paper I 
for details). We derived the enclosed Md yn ("IS <\ 
function of projected radius from position- velocity 
maps of the millimeter CO J = 1 — > 0, [Ne II] 
12.8 ^m and [S III] 9069A line emission (Shen & 
Lo 1995; Achtermann & Lacy 1995; McKeith et al. 
1993) combined with the result of stellar velocity 
dispersion measurements at the nucleus based on 
the CO 2.29 /um bandhead by Gaffney, Lester, & 
Telesco (1993). We estimated M gas , accounting for 
molecular and ionized gas, from the CO J = 1 — > 
map of Shen & Lo (1995) using the CO intensity 
to H 2 column density conversion factor derived by 
Wild et al. (1992) from multi-line radiative trans- 
fer analysis and from the Lyman continuum lumi- 
nosities. The M* for the nucleus and the starburst 
core are 7.9 x 10 7 and 6.1 x 10 8 M©, respectively. 
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We investigated the low-mass IMF (shape and 
cutoff) by comparing the mass of stars predicted 
by the model bursts M* tb with the M* derived 
for each region. In this comparison, allowance has 
to be made for a possible stellar population pre- 
dating the starburst, which can contribute to the 
measured stellar mass, so that the computed M* th 
does not exceed a plausible fraction of M*. To es- 
timate this fraction, we followed the argument of 
McLeod ct al. (1993). Tidal interaction with M 81 
presumably induced the starburst in M 82; numer- 
ical simulations indicate that the perturbations re- 
sulting from such encounters can efficiently drive 
large amounts of gas towards the nuclear regions 
and that before the gas mass equals the mass in 
preexisting stars, self-gravitation in the gas should 
trigger a starburst with short timescale compared 
to the gas settling time (e.g., Hernquist 1989; Mi- 
hos & Hernquist 1994, 1996). Thus, M* tb unlikely 
exceeds half of Md yn in M82. In our analysis, we 
adopted a conservative upper limit for M* tb of 50% 
of the total stellar mass M*. We did not account 
for the gas mass returned to the ISM via stellar 
winds and supernova explosions, or expelled out 
of the galaxy in the starburst wind. Since mas- 
sive stars contribute most to the mass return for 
the burst ages < 20 Myr relevant here but do not 
dominate M* tb , these effects can be neglected to 
first order but will be further addressed below. 

Our two-burst models for the nucleus and the 
starburst core (Tables 3 and 4) use only 9% 
and 33%, respectively, of the stellar masses for 
a Salpeter IMF slope between m\ ow — 1 M© and 
Wup = 100 M©, but this choice of m\ ow was arbi- 
trary. Figure 9 shows the burst masses required 
to reproduce the observed luminosities (Lk and 
-^Lyc) computed for m\ ow in the range 0.1 — 5 M©, 
with a Salpeter slope and for the same burst ages 
and decay timescales. For the nucleus, the bursts 
never use up more than 30% of the measured M*. 
On the other hand, the burst mass for the star- 
burst core reaches 50% of M* at m\ ov , = 0.4 M 
and equals M* at mi ow w 0.1 M©. For the mass in- 
volved in the starburst not to exceed our adopted 
limit would require a sharp cutoff at 0.4 M© for 
a Salpeter IMF. Alternatively, allowing the IMF 
to extend down to 0.1 M©, a flattening of the 
IMF at low masses would be required to reduce 
the increase in M* th with decreasing mi ow so that 
M* b < 50% of M*. With only one constraint 



available, our analysis does not allow us to break 
the degeneracy between the cutoff and the shape 
of the low-mass IMF. 

4-6.2. Factors Influencing the Low-mass IMF 
Determination 

The M* tb predicted by the models is deter- 
mined by the absolute luminosities. This makes 
inferences on the low-mass IMF very sensitive to 
the extinction corrections applied, to the star for- 
mation history, and to the shape of the IMF at 
higher masses. The first point has been exten- 
sively discussed by McLeod et al. (1993), Rieke 
ct al. (1993), and Satyapal et al. (1997). In that 
respect, our intrinsic luminosities correspond to 
the higher ranges reported in the literature and 
therefore impose more stringent constraints. In 
particular, the Lk used by Satyapal et al. (1997) 
and the lower limit considered by Rieke et al. 
(1993) are both smaller than our value (5 x 10 8 
and 7.8 x 10 s L© compared to 1.3 x 10 9 L© for a 
A = 1.9 — 2.5 fim K bandpass and the photomet- 
ric system of Wamsteker 1981 we adopted). With 
our two-burst models, they allow m\ ow down to 
0.2 M© or less, as shown in the bottom panel of 
Figure 9 (the mass constraints are similar between 
all three studies). 

The star formation history can also signifi- 
cantly alter the conclusions about the low-mass 
IMF since longer or older bursts tend to use up 
more mass. In this regard, however, our conclu- 
sions above for a Salpeter IMF differ from those of 
Satyapal et al. (1997) mainly because of the dif- 
ferent Lk adopted. In terms of M* th , their single 
burst model for the starburst core, with an age 
of 10 Myr and a short duration, is comparable to 
our two-burst model in which the old burst with 
an age of 9 Myr contributes 75% of the total burst 
mass for this region. Comparison with the models 
of Rieke et al. (1993) is less straightforward; our 
two-burst model is very similar to theirs with re- 
spect to burst ages and durations but the adopted 
IMFs are very different. 

For the young ages of interest here, the various 
luminosities are dominated by high-mass stars, so 
that the relative slope of the IMF between dif- 
ferent mass ranges is particularly critical in con- 
straining the low-mass IMF. For instance, IMFs 
which steepen at high masses are less efficient 
at converting the mass into luminosity. In order 
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to reproduce the luminosities, they require larger 
M* tb . To illustrate this, we considered three alter- 
native IMFs: the Miller & Scalo (1979) IMF for 
the solar neighbourhood, the IMF proposed by 
Rieke et al. (1993) for M82 (their "IMF 8"), and 
the IMF of Eisenhauer et al. (1998) for NGC 3603. 
Eisenhauer et al. derived the NGC 3603 IMF 
based on near-IR adaptive optics observations 
down to 1 M Q , their detection limit. For our 
purposes, we allowed extensions below 1 M us- 
ing the slope determined in the 1 — 15 M© range. 
Table 5 gives the parameters of the IMFs (see also 
Figure 8) and Figure 9 shows the variations of the 
burst masses versus m\ ow computed with each of 
them. 

For the range 0.1 — 100 M , the mass fraction 
locked in stars less massive than 1 M is signif- 
icantly smaller for the Miller & Scalo IMF com- 
pared to the Salpeter IMF (44% and 61%), while 
the corresponding fractions for stars more massive 
than 10 M are nearly identical (8% and 12%, re- 
spectively). However, properties such as Lk, £boi, 
and LLyc are steeply increasing functions of the 
stellar mass and the Miller & Scalo IMF is much 
steeper at the high-mass end than the Salpeter 
IMF. Consequently, normalizing to the same lu- 
minosities, the Miller & Scalo IMF predicts sub- 
stantially more mass than the Salpeter IMF, with 
M* tb larger by up to a factor of two in our mod- 
els. The Rieke et al. (1993) IMF is very similar to 
the Miller & Scalo IMF but with the first inflec- 
tion point moved up from 1 M , giving particular 
emphasis on the intermediate-mass range. Of the 
various IMFs they explored, such modifications 
to solar-neighbourhood type IMFs were the most 
successful at reproducing their set of constraints. 
Finally, compared to the Rieke et al. (1993) IMF, 
that of Eisenhauer et al. (1998) emphasizes less the 
intermediate-mass range but has a larger propor- 
tion of stars above 10 M Q (40% versus 23% of the 
total mass between 1 and 100 M Q ), implying much 
higher mass-to-luminosity conversion efficiency. 

Figure 9 shows that for the central 35 pc of 
M 82, the low-mass IMF remains unconstrained for 
any of the three alternative IMFs: the computed 
M* tb never exceeds 35% of the measured M* with 
toi ow down to 0.1 M . The case of the starburst 
core is again more restrictive. With the Miller & 
Scalo (1979) IMF, more than 50% of the stellar 
mass is involved in the bursts for m\ ow < 2 M . If 



we allow for stars down to 0.1 M to form, we con- 
clude as Rieke et al. (1993) that a displacement of 
the first inflection point to higher masses is then 
needed to effectively reduce M* th . With our data 
and models, though, M* th still exceeds half of M* 
at mi ow < 1 M if we adopt the Rieke et al. IMF. 
The Eisenhauer et al. (1998) IMF allows us to ex- 
tend the IMF with the same a — 1.73 power-law 
index from 15 M Q down to 0.1 M without us- 
ing up more than 30% of the stellar mass in the 
bursts. Interestingly, if we assume both the IMF 
and the if-band luminosity of Rieke et al. (1993), 
the M* tb remains well below 50% of M* over the 
entire range of m\ ov , (dashed line in the bottom 
panel of Figure 9). 

Another consideration worth pointing out is the 
mass loss via massive stars winds and supernova 
explosions, which is important already within a 
few million years of the onset of star formation 
(see, e.g., Lcithcrcr et al. 1999). In particular, 
the mass returned to the ISM could have been 
entrained out of M82 in the starburst wind, in 
which case our M* constraint would be missing 
some fraction or, alternatively, the predicted M* th 
curves of Figure 9 would decrease accordingly. An 
estimate of the effect is difficult because of the 
unknown amount of returned mass actually driven 
away in the starburst wind. As an illustration, we 
may consider the extreme case in which massive 
stars > 10 M lose all their mass within w 20 Myr, 
and that all ejectae are carried out of M82 by 
the starburst wind. For a Salpeter IMF extending 
from 1 to 100 M , the range > 10 M accounts 
for 31% of the total mass so that at most ss 30% 
of M* th would be lost and expelled. This fraction 
varies with mi ow , up to w 50% at 3 M and down 
to « 10% at 0.1 M , and more generally depends 
on all IMF parameters. The overall impact is to 
relax the requirements on the low-mass IMF but 
the effect is not likely to be more important than 
the other factors discussed above. 

4-6.3. Additional Remarks 

In view of the uncertainties involved and the 
lack of strong constraints on the shape of the IMF, 
the low-mass IMF in M 82 remains elusive and 
only general statements can be made. If stars 
formed down to 0.1 M , our data and models 
may indicate a substantial flattening of the star- 
burst IMF below a few M , assuming a Salpeter 
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slope at higher masses. This requirement would 
be alleviated for a flatter intermediate- and high- 
mass IMF or for less conservative limits on M* th . 
A bias against the formation of low-mass stars 
has been suggested for other starburst galaxies 
as well (e.g., Augarde & Lequeux 1985; Wright 
et al. 1988; Olofsson 1989; Nakagawa et al. 1989; 
Prestwich, Joseph, & Wright 1994). Engelbracht 
et al. (1996, 1998) concluded for NGC6946 and 
NGC 253 that the IMF proposed by Rieke et al. 
(1993) for M82 is more compatible with the ob- 
servations than a solar-neighbourhood IMF, while 
Alonso-Herrero et al. (2001) found that the IMF 
in NGC 1614 may be even more biased towards 
high-mass stars. 

Kinematic mass-to-light ratios have been ob- 
tained for a few individual luminous stellar clus- 
ters in M82, from velocity dispersion measure- 
ments combined with high resolution Hubble Space 
Telescope (HST) imaging. Using optical data ob- 
tained for the super star cluster "M 82-F" « 450 pc 
southwest of the nucleus and an age estimate of 
60 Myr, Smith & Gallagher (2001) concluded that 
the cluster IMF is top-heavy and requires a low- 
mass cutoff at 2 — 3 M for an a = 2.3 slope. 
McCrady, Gilbert, & Graham (2003) studied two 
younger (~ 10 7 yr) clusters roughly 180 pc south- 
west of the nucleus based on near-IR observations. 
For an IMF with a = 1.3 at 0.1 - 0.5 M Q and 
a = 2.3 at 0.5 - 100 M Q (based on Kroupa 2001), 
their data support a relative deficiency at < 1 M 
for one cluster but allow mi ow down to 0.1 M Q for 
the other. In their revision of the Smith & Gal- 
lagher (2001) analysis of M82-F, using the same 
mass but combined with the near-IR luminosity 
and a lower age of 40 Myr, they confirmed the ev- 
idence for a top-heavy IMF although with a looser 
requirement of m\ ow ss 1 M Q . Similar investiga- 
tions of young super star clusters in other systems, 
including NGC 1569A and NGC 1701-1 (Sternberg 
1998), and five clusters in NGC 4038/4039 (Men- 
gel et al. 2002), yielded mixed results as well. This 
seems to indicate a range of possible IMFs in su- 
per star clusters which may be related to their 
birth environment or to mass segregation (see, 
e.g., Mengel et al. 2002; McCrady et al. 2003). 
The results of the M 82 clusters may provide sup- 
port in favour of an IMF flattening at low masses 
but, as pointed out by Smith & Gallagher (2001), 
it is unclear whether conclusions about luminous 



star clusters, let alone of a handful only, can be 
extended to the overall pattern of star formation 
in the starburst. 

Direct censuses of the low-mass star popula- 
tions in Galactic and near-extragalactic high-mass 
star-forming regions are still scarce and difficult 
owing to the faintness of such stars and to impor- 
tant crowding effects. Recent results from very 
high sensitivity and angular resolution observa- 
tions suggest that mere truncations of the IMF 
near 1 M Q are not plausible, but the shape of the 
low-mass IMF remains poorly constrained. In par- 
ticular, Brandl et al. (1999) found that the cen- 
tral parsec of NGC 3603 is well populated by pre- 
main sequence stars down to 0.1 M Q , based on 
data from the Very Large Telescope. The near-IR 
adaptive optics observations by Eisenhauer et al. 
(1998) show no evidence for a significant flatten- 
ing of the IMF down to their detection limit of 
1 M . In contrast, from deep HST images of the 
R136 cluster in 30 Doradus, Sirianni et al. (1998, 
2000) could probe the 1.35 — 6.5 M Q range and 
reported a flattening of the IMF below 2 M , a 
result however disputed by Selman et al. (1999). 
Given the shape of the solar-neighbourhood IMF, 
with possible turnover near 0.5 M Q (e.g., Scalo 
1986, 1998; Rana 1987; Kroupa 2001), and indica- 
tions of a flattening towards lower masses in some 
local templates of starburst regions and super star 
clusters, our inferences for M 82 may not need be 
interpreted in terms of an "abnormal" deficiency 
of the IMF in low-mass stars. 

5. SPATIALLY DETAILED MODELING 
OF M82 

The modeling of selected regions in the previ- 
ous section provides the general characteristics of 
starburst activity in the central regions of M82, 
but only an incomplete picture of its evolution. 
In particular, our choice of the small-scale regions 
as those with brightest continuum or line emis- 
sion, and with deepest or shallowest CO band- 
heads likely introduced a bias towards preferen- 
tial ages. In order to obtain a more complete 
picture, we modeled individual regions throughout 
the entire 3D field of view. These have a range in 
their properties suggesting a range in evolutionary 
states and burst strengths. Although along any 
line of sight the integrated properties will always 
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be dominated by the most luminous populations, 
a detailed modeling on small spatial scales reduces 
the bias towards particular ages. 

Given our results of § 4, we assumed that a 
Salpeter IMF with high upper mass limits and 
burst decay timescales < 5 Myr are appropriate 
for individual pixels. We also assumed that the 
near-IR properties are dominated everywhere by 
red supergiants. We emphasize again that the 
analysis of the M*/Lk ratio of paper I strongly 
supports this hypothesis for the regions where the 
CO bandheads alone do not allow us to constrain 
the luminosity class of the evolved stars. 

5.1. Spatial Distribution of the Burst Ages 
and Strengths 

Figure 10 compares the data with single-burst 
models for m up = 100 M , and t sc = 1 and 5 Myr, 
in the diagrams of neon, iboi/^Lyc, and Lk/L^c 
ratios versus W2.29- The size of the data points 
is proportional to the intrinsic Lk- We did not 
account for the contribution of the cool evolved 
stars to the bolomctric luminosity of individual 
pixels, so that the values of L^/L\, yc are plotted 
in the iboi/^Lyc diagram. 

A large fraction of the data points in all three 
diagrams do not fall on the model curves and ex- 
hibit the typical behaviour encountered for the 
selected regions: the neon and Lg§/L Lyc ra- 
tios imply significantly younger ages than W2.29, 
with Lk/L^c corresponding to intermediate ages. 
In the Lk/Ll yc versus W2.29 diagram, the data 
points are distributed along a path that is remark- 
ably parallel to the model curve for t sc — 1 Myr, 
but displaced towards lower Lk / ^Lyc by about an 
order of magnitude. Although the properties of 
some of the individual pixels could be reconciled 
with a single burst if the timescale is increased to 
values approaching 5 Myr, we have pursued with 
our hypothesis of sequential bursts. The data can 
be reproduced with two short successive bursts, 
with the young bursts producing about ten times 
more ionizing luminosity than the older bursts 
dominating the near-IR luminosity. 

For each pixel, we therefore fitted a two-burst 
model to the constraints, each with t sc = 1 Myr, 
in the same manner as for the selected regions. 
We assigned an age of 1 Myr for Lj^ / ^Lyc ratios 
lower than the minimum Lboi/iLyc predicted by 



the models (this is the case for a few pixels only). 
For the old bursts, we adopted the youngest so- 
lutions inferred from W2.29 and assigned the age 
of maximum W2.29 to the few pixels with EWs ex- 
ceeding this limit. We estimated the formal uncer- 
tainties on the burst ages and initial star formation 
rates i?o by varying in turn the observational con- 
straints in the ranges given by the measurements 
uncertainties. 

Figure 11 presents the results. The ages for the 
young bursts range from 3.5 to 6.6 Myr, with an 
average of 4.8 Myr and dispersion of <r = 0.5 Myr. 
The ages for the old bursts range from 7.2 to 
12.6 Myr (the upper limit from the maximum 
W2.29), with an average of 8.8 Myr and disper- 
sion of (7 = 0.9 Myr. Figure 12 shows the average 
burst age and total burst strength as a function 
of projected distance from the nucleus in a slit 
3"-wide along the galactic plane of M82, as indi- 
cated in Figure 11. The typical uncertainties are 
±0.5 — 1 Myr for the ages and about 30% for the 
strengths. In view of these uncertainties, we con- 
clude that the maps and profiles indicate roughly 
constant ages among the young bursts and among 
the old bursts, as well as nearly uniform strengths 
for the young bursts. Spatial variations in strength 
among the old bursts are more significant with, in 
particular, a clear peak around the nucleus con- 
sistent with our findings from the modeling of se- 
lected regions (§ 4.4). 

Satyapal et al. (1997) investigated the spatial 
variations in ages from a dozen compact X-band 
emission sources, most of which are located within 
the regions mapped with 3D. They inferred burst 
ages in the range 4—10 Myr, with a trend of 
younger ages at larger projected radii. Given 
that they modeled the sources as single bursts 
and used the strength of the CO bandheads long- 
wards of 2.29 /im together with the EW of the Br7 
emission line at 2.17^m (inversely proportional to 
Lk/Llyc), their results are consistent with ours. 

5.2. The Global Star Formation History 

From our spatially detailed modeling, we recon- 
structed the global star formation history within 
the 3D field of view by integrating the initial star 
formation rate Rq as a function of burst age th 
over all individual pixels. For sufficient sampling, 
we chose age bins increasing logarithmically by 
0.05 dex. As long as the timescales assumed for all 
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individual bursts are similar and relatively short 
compared to the overall star formation history of 
interest, the integrated i?o(£b) curve gives a mea- 
sure of the "instantaneous" star formation rate at 
different times. Figure 13 shows the derived inte- 
grated starburst history together with the surface 
density of star formation rate, computed by divid- 
ing the integrated i?o(^b) by the total area of the 
pixels contributing to each age bin. The results 
for the individual pixels are reproduced as well. 

The integrated star formation rate exhibits two 
conspicuous peaks near 5 and 9 Myr, correspond- 
ing to distinct global starburst episodes. The first 
episode was globally 2.5 times stronger than the 
most recent one in terms of mass of stars formed. 
The star formation rate surface density reflects the 
nearly uniform strengths for the young bursts and 
the larger variations for the old bursts seen in Fig- 
ures 11 and 12. We note that modeling the indi- 
vidual pixels with two short successive bursts in- 
troduces a bias towards separated burst events in 
the integrated star formation history. Again, such 
a scenario is however plausible in view of the spa- 
tial anticorrelation between the tracers of OB stars 
and those of cool evolved stars (§ 4.1). 

The integrated star formation rate in the 3D 
field of view is well reproduced by two Gaussians 
in Rq versus age centered at 4.7 and 8.9 Myr, 
with amplitudes of 6.3 and 18.5 M Q yr _1 and full- 
widths at half-maximum of 1.1 and 1.7 Myr, re- 
spectively. These Gaussians are shown in Fig- 
ure 13 along with the two-burst model fit to the 
integrated properties of the 3D field of view from 
§ 4. The exponentially decaying functions pro- 
vide reasonable though crude approximations to 
the detailed star formation history. In view of our 
results, it may not be surprising that Rieke et al. 
(1993) found that a double Gaussian profile opti- 
mized their model fits for the entire starburst core 
(with one of their best models reaching peak star 
formation activity near 5 and 13 Myr). 

6. THE NATURE OF STARBURST AC- 
TIVITY IN M 82 

6.1. Star Formation Process 

Observations of M82 reveal important small- 
scale structure tracing individual burst sites on 
10-pc or even parsec scales (see paper I). In par- 
ticular, the ionized and molecular gas as well as 



the IR-emitting dust exhibit clumpy morphologies 
on scales at least as small as 15 pc (e.g., Telesco 
& Gezari 1992; Larkin et al. 1994; Achtermann & 
Lacy 1995; Shcn & Lo 1995). Compact continuum 
sources tracing clusters of red supergiants near 
their maximum luminosity are detected at near-IR 
wavelengths (e.g., Satyapal et al. 1997; McCrady 
et al. 2003). Optical imaging with HST has re- 
solved over a hundred young super star clusters, 
some of which may lie in directions of lower ex- 
tinction and belong to the inner starburst regions 
(O'Connell et al. 1995). 

Our models show that the properties of individ- 
ual burst sites are consistent with the formation of 
very massive stars and with short burst timescales 
of a few million years or less (see also Satyapal et 
al. 1997, for selected clusters). The typical sizes 
and the star formation parameters inferred make 
these burst sites comparable to Galactic and near- 
extragalactic massive star-forming regions such as 
NGC3603 and the R136 cluster in 30 Doradus. 

Interestingly, the smoother and low-surface 
brightness emission regions make a substantial 
contribution to the integrated properties, suggest- 
ing that the prominent sub-structure may in fact 
trace the largest and most luminous burst sites. 
The diffuse ionized gas emission (notably in H 
recombination lines) and thermal dust emission 
represent about 50% of the respective total emis- 
sion from the starburst core. The smoother and 
fainter stellar if-band continuum emission con- 
tributes at least as large a fraction to the total 
if-band luminosity. These components corre- 
late roughly with the overall distribution of the 
brighter and more compact sources, and partly 
break up into smaller-scale structure at higher 
angular resolution (see, e.g., the HST NICMOS 
maps of Alonso-Hcrrcro et al. 2003). Whether the 
diffuse emission in M 82, and more generally in 
starburst galaxies, traces large numbers of unre- 
solved fainter compact sources or truly extended 
sources is still an open question, with implications 
for the importance of cluster formation in star- 
bursts and subsequent dissolution in the stellar 
field (e.g., Ho & Fillipenko 1996; Meurer et al. 
1995; Tremonti et al. 2001). 



14 



6.2. Evolution of Starburst Activity 

6.2.1. Key Morphological Features 

For a plausible evolutionary scenario of star- 
burst activity in M82, our results must be inter- 
preted together with the key morphological fea- 
tures: 

• Large-scale tails and bridges of material ema- 
nating from M 82 and connecting with its massive 
neighbour M 81 located about 36 kpc in projection 
indicate gravitational interaction between the two 
galaxies (e.g., Yun, Ho, & Lo 1993, 1994). 

• The morphology of the -fT-band emission, dom- 
inated by red supergiants, is highly suggestive of 
a nearly edge-on disk, peaks strongly towards the 
nucleus, and provides the strongest evidence for 
a ~ 1 kpc-long stellar bar (e.g., Telesco et al. 
1991; McLeod et al. 1993; Larkin et al. 1994). The 
brightest and most compact sources are mainly lo- 
cated in an "inner disk" within « 150 pc of the 
nucleus and at the ends of the bar (e.g., Satyapal 
et al. 1997). 

• The most prominent H II regions are concen- 
trated in a rotating ring-like structure of radius 
« 85 pc and, outside this ring, along ridges pre- 
sumably on the leading edge of the rotating stellar 
bar (e.g., Larkin et al. 1994; Achtermann & Lacy 
1995). The "ionized ring" is located just inside 
of the main concentrations of molecular gas in a 
toroid or tightly-wound spiral arms (e.g., Shen & 
Lo 1995). Our Br7 source B2 coincides in pro- 
jection with the western edge of the ionized ring 
while Bl is located farther out. 

• The important population of young supernova 
remnants (SNRs) detected at radio wavelengths 
is quite uniformly distributed along the galactic 
plane over « 600 pc (e.g., Kronberg, Biermann, 
& Schwab 1985; Muxlow et al. 1994). Their spa- 
tial distribution together with the ages derived 
for each starburst episode and the IMF weight- 
ing are more consistent with their being associated 
primarily with the populations of red supergiants 
throughout the disk. 

• A bipolar outflow along the minor axis of M 82 
is traced out to at least 5 kpc by X-ray and 
optical observations (e.g., Bregman, Schulman, 
& Tomisaka 1995; Shopbell & Bland-Hawthorn 
1998). The [Fe II] 1.644 ^m emission exhibits a 
prominent arc-like structure to the south and cen- 
tered near the nucleus. It may trace shock-ionized 



iron-enriched gas where the outflowing wind inter- 
acts with interstellar disk material (e.g., Green- 
house et al. 1997; Alonso-Herrero et al. 2003). 
• There is an apparent lack of gas and dust within 
the central few tens of parsecs of M82 (e.g., Tele- 
sco & Gezari 1992; Larkin et al. 1994; Shen & Lo 
1995; Achtermann & Lacy 1995; Seaquist et al. 
1996). 

Given the relative distributions of the gaseous 
and stellar components and since M82 is viewed 
at an inclination angle of « 80°, the two sequen- 
tial bursts inferred locally and globally within the 
central starburst regions could be understood in 
terms of projection effects. In particular, the inner 
few tens of parsecs harbored the most intense star- 
burst activity about 10 Myr ago and while more 
recent bursts appear superposed, they took place 
predominantly in circumnuclear regions at larger 
radii. The geometrical picture of M82 further em- 
phasizes the particular intensity of the old nuclear 
burst revealed in Figures 11 and 12, with the red 
supergiants within the central few tens of parsecs 
in projection being physically close to the nucleus 
along the line of sight as well, while the OB stars 
— and the red supergiants elsewhere — possibly 
more extended along the line of sight. While alter- 
natives could be advanced, we chose, in line with 
our two-burst hypothesis, to focus on the above 
interpretation as an illustrative yet possible case. 

6.2.2. Triggering Mechanism 

The primary triggering mechanism for star- 
burst activity in M82 is generally attributed to 
the M82 - M81 tidal interaction ~ 10 8 yr ago 
(e.g., Gottesman & Weliachew 1977; O'Connell & 
Mangano 1978; Lo et al. 1987; Yun et al. 1993, 
1994). In this scenario, the ISM in M82 experi- 
enced strong large-scale torques and loss of angu- 
lar momentum as it was transported towards the 
dynamical center of the galaxy, in accordance with 
numerical simulations (e.g., Sundelius et al. 1987; 
Noguchi 1987, 1988; Mihos & Hernquist 1996). 

The increased cloud-cloud collision rate in the 
disk and the large amounts of material accumu- 
lated and compressed in the innermost regions 
could have led to the first starburst episode, char- 
acterized by very intense star formation at the 
nucleus and lower-level activity elsewhere in the 
disk. The nuclear burst and subsequent high rate 
of supernova explosions consumed rapidly the gas 
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supply and expelled the remaining gas via the star- 
burst wind, thereby creating the central cavity in 
the ISM and preventing further star formation. 
The absence of a large concentration of radio SNRs 
in the nuclear vicinity is consistent with the lack 
of gas and dust. 

6.2.3. Subsequent Evolution 

Numerical simulations show that bars in galac- 
tic disks can be induced by galactic interactions 
and are very effective at driving material towards 
the central regions of galaxies (e.g., Combes & 
Gerin 1985; Noguchi 1987, 1988; Shlosman et al. 
1989; Athanassoula 1992). Ring- or spiral-like 
dynamical resonances may develop under the ac- 
tion of such non-axisymmetric perturbations. In 
the presence of inner Lindblad resonances (ILRs), 
the radial inflow of material is stopped before it 
reaches the nucleus of the galaxy and accumulates 
in a circumnuclcar ring. Star formation is trig- 
gered by shocks in the ring and along ridges lead- 
ing the bar, and may be particularly enhanced at 
their intersections. 

The second starburst episode in our models of 
M82 could thus be attributed to the presence of 
the bar and accompanying ILRs indicated by the 
global gas distributions (see also Lo et al. 1987; 
Telesco et al. 1991; Achtcrmann & Lacy 1995; 
Neininger et al. 1998). The ionized ridges pre- 
sumably trace star formation sites on the leading 
side of the bar. Enhanced activity where the ISM 
streaming along the bar meets with the ionized 
ring may be hinted at by the slight increase in 
strength of the young bursts near B2 suggested by 
Figures 11 and 12 (see also Achtcrmann & Lacy 
1995). The bar may have played a role during the 
first starburst episode by channeling the nuclear 
inflow before the present dynamical resonances ap- 
peared. 

The apparent outward progression of starburst 
activity suggested by the sequence Nuclear super- 
giants — » Ionized ring — > Molecular ring may re- 
flect the temporal development of the bar and re- 
distribution of the mass as suggested by Telesco et 
al. (1991). At the radius of the ionized ring, the 
rotational velocity is about 120 kms -1 (Achter- 
mann & Lacy 1995; Shen & Lo 1995). This implies 
an orbital period of 4 — 5 Myr, i.e. a dynamical 
timescale similar to the time elapsed between the 
peak star formation activity during the two recent 



starburst episodes. 

An alternative explanation could be self-induced 
propagation radially outward in the disk, as a 
consequence of the powerful nuclear burst and ex- 
panding shock wave generated by the massive stars 
winds and supernova explosions (e.g., Carlstrom 
& Kronberg 1991; Achtcrmann & Lacy 1995). 
Without dynamical resonances, the ISM may be 
expected to follow a smooth, radially decreasing 
density profile. Star formation triggered by an ex- 
panding "superbubblc" in such a medium seems 
difficult to reconcile with pronounced and narrow 
peaks in the global star formation rate. The two 
peaks could perhaps reflect a sudden increase of 
mechanical energy release by the older nuclear 
burst after about 5 Myr but the contribution by 
massive star winds is already very substantial at 
ages < 5 Myr before the first supernovae explode 
(e.g., Leitherer et al. 1999). Star formation ac- 
tivity along the ridges beyond the ionized ring is 
also difficult to explain in the self-induced outward 
propagation scenario. 

We favoured bar-driven evolution as dominant 
mechanism because the stellar and gaseous compo- 
nents exhibit morphological features strongly sup- 
porting the existence of a bar and its dynamical 
effects. Bar-driven evolution has been proposed 
for other starburst systems as well as more quies- 
cent spiral galaxies in which circumnuclear rings 
of enhanced gas density and star formation ac- 
tivity coexist with stellar or gaseous bars (e.g., 
Telesco & Decher 1988; Telesco, Dressel, & Wols- 
tencroft 1993; Knapen et al. 1995; Boker, Forster 
Schreiber, & Genzel 1997). 

6.24. Global Picture 

Investigations of regions outside of the central 
active starburst core outline further key elements 
for a global scenario. The integrated optical spec- 
trum of M 82 exhibits strong Balmer line transi- 
tions in absorption typical of A-type stars, indica- 
tive of an intermediate-age (0.1 — 1 Gyr) "post- 
starburst" population (Kennicutt 1992). de Grijs 
et al. (2000, 2001) used HST optical and near-IR 
data to study post-starburst regions 0.5 — 1 kpc 
northeast of the nucleus, identified as such from 
strong Balmer absorption lines and discontinuity 
in their spectra (O'Connell & Mangano 1978; Mar- 
cum & O'Connell 1996). The results of de Grijs 
and coworkers based on evolutionary synthesis 
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modeling of the photometric properties of compact 
luminous clusters and bright giant stars as well 
as their detection of candidate Ha SNRs imply 
that enhanced star formation occurred ~ 600 Myr 
ago at radii of ~ 1 kpc and up to ~ 30 Myr ago 
near - 500 pc. Gallagher & Smith (1999) mod- 
eled high-quality optical spectra of two super star 
clusters about 450 pc southwest of the nucleus and 
derived ages of <~ 60 Myr. 

In the framework of our models and interpreta- 
tion, together with the spatial distribution of the 
radio SNRs and the M* /Lk ratio, the first star- 
burst episode ~ 10 Myr ago was not strictly con- 
fined to the nuclear vicinity, but occurred through- 
out the central 500 pc as well. Along with the ev- 
idence from the post-starburst regions, this would 
support an "outside-inside-out" propagation of 
star burst activity in M82. The scenario we pro- 
pose is summarized here and sketched in Fig- 
ure 14. 

• Following the gravitational interaction be- 
tween M82 and M81, ~ 10 8 yr ago, the ISM in 
M 82 experienced strong large-scale torques, loss of 
angular momentum, and important infall towards 
the nuclear regions, leading to enhanced star for- 
mation activity in the central kiloparsec. The stel- 
lar bar induced by the interaction possibly played 
a role in channeling the inflow. In the absence of 
dynamical resonances, the infalling material can 
reach the nucleus. 

• Within the central 500 pc of M 82, a first star- 
burst episode took place 8—15 Myr ago and was 
rapidly exhausted. The inner few tens of parsecs 
at the nucleus hosted the most intense star forma- 
tion activity. 

• A subsequent starburst episode was trig- 
gered predominantly by bar-induced dynamical 
resonances. It occurred 4 — 6 Myr ago, mainly 
in a circumnuclear ring and along the stellar bar, 
and also decayed rapidly. 

• A supernova-driven starburst wind originat- 
ing in the center of M82 has broken out of the 
galactic plane, the dramatic aftermath of the pow- 
erful nuclear burst. The outflow component in the 
disk may have played some role in triggering the 
second starburst episode. 



6.3. Timescales and Feedback Effects 

The global durations of ^ 10 6 yr derived for 
each of the recent starburst episodes in M 82 seem 
surprising for spatial scales up to at least 500 pc 
and in view of the longer estimates ~ 10 7 — 10 8 yr 
often quoted for starburst systems (e.g., Thron- 
son & Telesco 1986; Heckman 1998). We consider 
below three timescale arguments in trying to gain 
more insight into the quenching mechanisms of the 
recent starbursts in M82. 

One argument relics on the comparison between 
the present star formation rate and the mass of 
the current gas reservoir, giving a "gas consump- 
tion timescale" T gas . In M82, the gas is mostly in 
molecular form, with a total M H2 = 1.8 x 10 8 M 
(e.g., Wild et al. 1992). From our two-burst model 
of the starburst core (Table 4) and with a Salpeter 
IMF, M* th « 2 x 10 8 M© of stars were formed 
in the past w 15 Myr for m\ ow — 1 M Q and 
w 5 x 10 8 Mq for m low = 0.1 M©. The latter 
value represents 80% of the measured stellar mass 
and, in view of the discussion of § 4.6, is here con- 
sidered as an upper limit. The average star for- 
mation rate in the recent history of M 82 is then 
R = 13 — 33 M Q yr _1 , depending on mi ow - As- 
suming 100% efficiency, starburst activity in M 82 
at the level observed in the past 15 Myr could be 
sustained for another r gas w 5 — 15 Myr. This 
provides a fair estimate for the minimum gas de- 
pletion timescale but has no connection with the 
physical processes responsible for the evolution of 
the starburst. 

Considerations based on dynamical timescales 
are more closely linked to the large-scale mech- 
anisms that may drive starburst activity. For a 
sample of 36 IR-selected starburst galaxies, Ken- 
nicutt (1998) inferred a median gas conversion ef- 
ficiency of 30% per dynamical timescale, taken 
as an average orbital period for all objects of 
10 8 yr; e gas = 100%/10 8 yr occurs but in the 
most extreme objects. From the rotation curve 
of M82 (e.g., Gotz et al. 1990; McKeith et al. 
1993; Achtermann & Lacy 1995; Shcn & Lo 1995; 
Neininger et al. 1998), we compute shorter dynam- 
ical timescales for regions inside the starburst core. 
For example, the orbital periods are about 5, 10, 
and 20 Myr at radii of 90, 225, and 500 pc cor- 
responding to the locations of the ionized ring, of 
the main concentrations of molecular gas, and of 
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the outer limit enclosing the starburst core and 
most of the gas. The dynamics allow then for 
burst timescales of approximately Td yn = 15 Myr, 
30 Myr, and 65 Myr for 30% efficiency per orbital 
period. Even with 100% efficiency, the timescales 
are barely short enough to explain the durations 
of each starburst episode in M82. 

Such dynamical considerations may provide 
more insight concerning the settling time of the 
gas and the development of dynamical resonances. 
The orbital time at the radius of the ionized ring 
is similar to the time separation between the two 
most recent bursts. At radii 1 — 3 kpc, the or- 
bital period is « 100 — 200 Myr, which is of the 
order of the age of the super star clusters in the 
post-starburst regions ~ 1 kpc northeast of the 
nucleus and to the epoch of last peri-passage be- 
tween M82 and M81. The scenario proposed in 
§ 6.2 for the triggering and propagation of the 
starburst episodes in M 82 is consistent with these 
timescales. 

None of the above arguments is however com- 
pletely satisfactory for explaining the very short 
global burst durations suggested by our mod- 
els. Most importantly, any feedback effect of 
massive stars via strong stellar winds and subse- 
quent supernova explosions is neglected. Qualita- 
tively, timescales of a few million years are sim- 
ilar to the lifetimes of stars with initial masses 
> 50 M . Although both episodes took place over 
~ 100 — 500 pc scales, each was triggered nearly 
simultaneously everywhere, presumably as a re- 
sult of brief gas compression events. Massive star 
winds and supernova explosions rapidly disrupt- 
ing the remaining ISM locally, preventing further 
star formation after a few million years, may have 
"conspired" to produce a rapid decay of the global 
starburst activity as well. 

We can make simple estimates for the "feed- 
back timescale" Tf ee d by comparing the cumulative 
mechanical energy -E mec h injected by supernovae 
(and neglecting the contribution by stellar winds) 
into the ISM with its gravitational binding energy 
-^bind, and asssuming that starburst activity stops 
when they exactly balance each other: 

-EmochTfocd = -^bind, (1) 

where -Emech is the rate of mechanical energy de- 
position. We assumed that all stars with initial 
masses > 8 M Q explode as supernovae with a 



typical energy of 10 51 erg, of which a fraction r\ 
is transferred as kinetic energy to the ISM. For 
simplicity, we considered a spherically symmetric 
distribution of gas and stars within the starburst 
core, so that -Ebind = (GMd yn M gas ) /1Z where G is 
the gravitational constant, Md yn and M gas are the 
dynamical and gas masses, and 1Z is the radius of 
the region. 

The rate of supernova explosions ^sn is the 
most direct tracer of i^mcch but the observed rates 
suffer from rather large uncertainties. We com- 
puted -Emcch using the inferred star formation rate 
R instead. The implicit assumption is that the 
burst is sufficiently evolved to have non-negligible 
i^sn and that the star formation rate is still sub- 
stantial. For a Salpeter IMF between 1 M and 
100 M , zajn [yr -1 ] w 0.02 R [Moyr" 1 ]. With 
these assumptions, Eq. 1 can be written: 

^£!~ 500 (L\ ( M ^ Yf—Y 1 ( - ) 
Myr U/UO 1O M ; \kpc J \M Q yv^J 

(2) 

where / is the gas mass fraction. For the starburst 
core ofM82, TZ = 225 pc, M dyn « 8xl0 8 M , and 
M gas w 1.8 x 10 8 M (see paper I). With the av- 
erage R m 13 M yr _1 derived for the appropriate 
IMF parameters, Eq. 2 yields rf co d — 0.27/ -1 Myr. 
The supernova efficiency is poorly constrained, 
but values in the range r) = 0.1 — 1 seem appro- 
priate for M82 (Chevalier & Clegg 1985). Con- 
sequently, Tf oc d — 0.2 — 2 Myr. For compari- 
son, with j^sn ~ 0.06 yr _1 for M82, we obtain 
Tfocd — 1 — 10 Myr depending on r\. As alternative 
tracer, we can use Lboi- In a steady-state, our star- 
burst models imply Lboi [L©] ~ 10 12 ^sn [y r_1 ], 
and thus R [M yr" 1 ] - 0.5 (iboi/10 10 L ). With 
£boi = 6.6 x 10 10 L for the starburst core, Eq. 2 
implies Tf oc d ~ 1 — 10 Myr as well. These estimates 
represent upper limits because we neglected mas- 
sive stars winds. 

The Tfecd is more consistent than T gas or Td yn 
with the short durations in our models for the re- 
cent starburst episodes in M82. Admittedly, our 
Tfocd estimates are also simplistic but they provide 
an additional perspective for the evolution of star- 
burst activity and a possible explanation for short 
global burst durations. Along with our spatially 
detailed models, the above timescales are consis- 
tent with the overall progression of starburst ac- 
tivity in M 82 being driven by large-scale dynam- 
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ical processes related to the tidal interaction with 
M 81 and the induced stellar bar, but suggests that 
once starburst activity was triggered, strong neg- 
ative feedback effects acted very rapidly to inhibit 
further star formation on all spatial scales. 

Recent studies of other starburst systems also 
provide evidence for short burst durations and 
for the episodic nature of starburst activity (e.g., 
Thornley et al. 2000; Boker et al. 1997; Vanzi et al. 
1998; Alonso-Hcrrero et al. 2001). Some numerical 
simulations of the dynamical evolution of barred 
galaxies which account for the effects of star for- 
mation indicate a recurrent burst behaviour with 
typical timescale ~ 10 7 yr, owing to the ability of 
massive stars to destroy the flow pattern when a 
sufficient number coexist in a given place (Heller 
& Shlosman 1994; Knapen et al. 1995). 

7. SUMMARY 

We have presented the results of detailed mod- 
cling of the central starburst regions of M82, 
based on new near-IR integral field spectroscopy 
and mid-IR spectroscopy complemented with ad- 
ditional results taken from the literature. We ap- 
plied starburst models optimized for observations 
in the range A = 1 — 45 /im to these data and con- 
strained the star formation parameters (cutoffs of 
the IMF, burst timescale) as well as the spatial and 
temporal evolution of starburst activity on scales 
as small as 25 pc. 

Our data and models are consistent with the 
formation of very massive stars (> 50 M Q ) and im- 
ply burst durations of at most a few million years 
for individual burst sites. With typical sizes on 
10-pc scales or less, individual burst sites in M 82 
are comparable to Galactic and near-extragalactic 
massive star-forming regions such as NGC3603 
and the R136 cluster in 30 Doradus. For plau- 
sible limits on the mass involved in the starburst, 
we find that the IMF must flatten at low masses 
assuming a Salpeter slope at high masses. Dynam- 
ical information for three individual young stellar 
clusters in M82 (Smith & Gallagher 2001; Mc- 
Crady et al. 2003) provide similar evidence, but 
clearly more clusters need to be investigated. In 
view of the large uncertainties involved and the 
scarcity of low-mass IMF determinations in local 
high-mass star-forming regions, this result is diffi- 
cult to interpret. 



Among possible alternatives, our model results 
and interpretation together with the spatial dis- 
tribution of the stellar and gaseous components in 
M 82 leads us to propose the following scenario. 
Globally, starburst activity in the central 500 pc 
of M82 occurred during two successive episodes 
each lasting a few million years, about 10 and 
5 Myr ago. The first episode took place through- 
out the disk and was particularly intense in the 
central few tens of parsecs. The second episode 
occurred predominantly in circumnuclear regions 
and along the stellar bar. We interpret this se- 
quence in a tidally-triggered, bar-driven evolution 
scenario, consistent with the evidence for gravita- 
tional interaction between M 82 and its neighbour 
M81, with the observed morphological signatures 
of the bar and associated dynamical resonances, 
and with the corresponding dynamical timescales. 

The episodic nature of starburst activity and 
short global burst decay timescales in M82 are 
particularly interesting results. Invoking a sim- 
ple argument based on the comparison of the cu- 
mulative mechanical energy injected into the ISM 
and its gravitational binding energy, we find that 
the collective effect of massive stars winds and su- 
pernova explosions can disrupt the ISM in M82 
on timescales of ~ 10 6 yr, providing a very ef- 
ficient quenching mechanism for each starburst 
episode. Following brief gas compression events 
on large-scales, starburst activity decayed very 
rapidly due to its own strong negative feedback ef- 
fects. Our scenario for M 82 outlines the interplay 
between the large-scale triggering mechanisms and 
the more local but important feedback effects in 
determining the evolution of starburst activity in 
this galaxy. 
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A. FURTHER DETAILS ON STARS 



A.l. Spectral Energy Distributions 

STARS computes the detailed SED intended for photoionization modeling of the nebular emission excited 
by the stars. The SEDs are taken from a hybrid grid generated as described by Thornley et al. (2000) 
from the LTE models of Kurucz (1992) for effective temperatures T c g < 19, 000 K and the non-LTE models 
of Pauldrach et al. (1998) for T e g > 25,000 K, with SEDs for intermediate temperatures obtained by 
interpolation. The Pauldrach et al. models represent much better the effects of line blocking and blanketing 
in the rapidly expanding atmospheres of hot stars and are generally harder than the Kurucz models. A 
detailed description of the models and implications on the computed SEDs is given by Pauldrach et al. 
(2001) and references therein. The basic Pauldrach et al. model SEDs we used in STARS include dwarfs 
and supergiants. For dwarfs, they cover 25000 — 60000 K in T e g at fixed surface gravity log(g) = 4 and 
mass-loss rates M of 10~ 8 - 5 x 10~ 6 Mgyr 1 . The ranges for supergiants are T eS = 25000 - 50000 K, 
log(g) = 2.75 - 3.8, and M of 2 x 10~ 6 - 1.5 x 10~ 5 M Q yr _1 . For a given T cff and L hol along a stellar track, 
STARS uses the corresponding evolutionary log(g) and assigns the SED model with closest T e g and log(g) 
from the interpolated hybrid grid. The SED model grid covers well the range in parameters of the Geneva 
stellar tracks during most of the hot stars lifetimes (for both normal and enhanced mass-loss rates sets). 

Various quantities predicted by STARS are derived using the empirical bolometric corrections of Schmidt- 
Kaler (1982) and the broad-band colours of Koornneef (1983). This is especially important for cool stars 
which have SEDs severely distorted by extensive molecular absorption and contaminated by thermal emission 
from a circumstellar dust shell that are not reproduced satisfactorily by black-body approximations or 
currently available model atmospheres. The impact of using empirical data is largest at near-IR wavelengths. 
For instance, the predicted if -band luminosities are up to 30% larger for the age range ~ 10 — 100 Myr 
(assuming t sc = 1 Myr). 

A. 2. Near-IR Stellar Absorption Features 

STARS predicts the EWs of near-IR stellar absorption features characteristic of late-type stars: the 
12 CO(6,3) and 12 CO(2,0) bandheads at 1.62 and 2.29 ^m, and the Si I, Na I, and Ca I features at 1.59, 
2.21, and 2.26 /im. We implemented in STARS the EWs versus T c ff relationships for luminosity classes V, 
III, and I below 5700 K obtained by polynomial fitting to the homogeneized stellar data compiled by Forster 
Schreiber (2000). The EWs assigned to hotter stars were identically zero since the corresponding features 
vanish at higher T c tf . The exception is the 1.59 /xm feature but above w 6000 K, the hydrogen Brl4 transition 
is responsible for the absorption feature which has then a different physical meaning. The predicted EWs 
include "dilution" effects due to featureless continuum emission from OB stars and from associated free-free 
and free-bound processes. The stellar data are given for spectral resolutions of R ~ 1600 and 2500 in the 
H and K band, respectively. For comparison with the models, EWs measured from lower resolution data 
should be corrected according to the relationships given by Forster Schreiber (2000). 

Recently, Origlia et al. (1999) and Origlia & Oliva (2000) questioned the reliability of synthesis modeling of 
the CO bandheads from a critical analysis of theoretical modeling and observational data. They emphasized 
the importance of accounting for AGB evolution up to the end of the thermally-pulsing phase. They also 
demonstrated that uncertainties in stellar evolutionary tracks for red supergiant and AGB phases remain 
large especially at sub-solar metallicities but affect much less the results for ages up to ~ 10 7 yr at near-solar 
metallicities. Since STARS accounts for the thermally-pulsing AGB (see below) and with the solar-metallicity 
tracks adopted in this work, these uncertainties are minimized in our predictions of the CO EWs. 

A.3. Implementation of the TP- AGB Phase 

The Geneva tracks (Schaller et al. 1992) employed in STARS follow the evolution of stars up to well defined 
evolutionary phases. Tracks for massive stars (M > 7 M Q ) extend up to the end of the core carbon-burning 
phase, just prior to their explosion as core-collapse supernovae. Tracks for low-mass stars (M < 2 M Q ) end 
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at the helium flash which occurs at ages > 10 9 yr. The subsequent evolution of low-mass stars is not relevant 
for the starburst ages (< 10 9 yr) we consider. For intermediate-mass stars (2 M < M < 7 M Q ), the Geneva 
tracks are terminated at the end of the early asymptotic giant branch (E-AGB), at ages in the range 5 x 10 7 
to 10 9 yr, and thus exclude the TP-AGB phase. TP-AGB stars can contribute significantly to the near-IR 
and total luminosities in starbursts, and so should be included in synthesis models (e.g., Chariot & Bruzual 
1991; Maraston 1998; Origlia & Oliva 2000; Mouhcine & Lancon 2002). Observationally, there is growing 
evidence that such intermediate-age populations account for a non-negligible fraction of the integrated light 
in the central regions or bulges of several nearby galaxies, including our own (e.g., Freedman 1992; Blum, 
Sellgrcn, & DcPoy 1996; McLean & Liu 1996; Lee 1996; Davidge et al. 1997; Martmez-Delgado & Aparicio 
1998). 

We extended the Geneva tracks for intermediate-mass stars using an approach suggested by Chariot & 
Bruzual (1991) which employs the analysis of the TP-AGB presented by Bedijn (1988). In this procedure, 
the TP-AGB is divided into a beginning "Mira phase" where the stars pulsate in the first overtone (or 
higher pulsation modes) while losing mass slowly, and a late "OH/IR" phase where the stars pulsate in the 
fundamental mode and lose mass rapidly. In the Mira phase, the stars increase in luminosity (and cool) 
and move up the AGB track appropriate for their mass. In the OH/IR phase, the luminosities do not vary 
much because of the rapid mass loss. Bedijn (1988) provides analysis and computations of the durations of 
each phase, including the maximum luminosities L max at the end of the OH/IR phase, as functions of the 
main-sequence masses between 2 and 7 M Q . Using this information, we added two evolutionary points in 
the Lboi versus T e g Hertzsprung-Russell diagram. The first point corresponds to the end of the Mira phase 
where (following Bedijn) we assume that Lboi = 0.8L max . The second point corresponds to the end of the 
OH/IR phase where Lboi = L max . We obtained the effective temperature at each point from the Lboi — Toff' 
relations for stars on the AGB. Our extensions are illustrated in Figure 15. As sanity check, we compared 
the Lboi predicted by our shortest burst models based on conventional synthesis and by applying the fuel 
consumption theorem (Renzini & Buzzoni 1986) to our library of tracks. We obtain similar values within 
10% on average up to 10 9 yr, with maximum differences of ±25% at < 3 x 10 8 yr increasing to ±40% at 
older ages where the numerical fluctuations due to the finite number of tracks become more important. 

In attributing near-IR properties to TP-AGB stars, we followed a simple approach which can be justified 
by the small dispersion in T c g of TP-AGB stars of various masses and the fact that the time spent in the 
OH/IR phase is at most 30% of the total TP-AGB lifetime according to Bedijn's models (see also Chariot & 
Bruzual 1991). We assigned to all TP-AGB stars the average bolometric correction and broad-band colours 
of the prototypical Mira variable u> Cet, derived from the data of Mendoza & Johnson (1965) and Zhou 
et al. (1984). For the EWs of near-IR absorption features, we adopted preliminary data obtained with the 
3D instrument for a sample of Miras and other TP-AGB stars (S- and N-type stars). The observed CO 
bandhead strengths are comparable to those of the latest normal M-type giants. Similar conclusions are 
reached by inspection of the spectra published, e.g., by Lancon et al. (1999) and Lancon & Wood (2000). 
We therefore attributed the EWs of normal giants with T off = 3000 K to all TP-AGB stars. 

As expected, the effects of including the TP-AGB phase are most important for Lboi and for the near-IR 
properties while they are negligible for optical and UV properties. Figure 16 shows the TP-AGB contribution 
in our synthetized Lbob Lk, and Detailed comparison with other models including the TP-AGB, 

notably by Chariot & Bruzual (1991), Maraston (1998), and Mouhcine & Lancon (2002), is difficult be- 
cause of the differences in synthesis techniques (conventional, isochrone, fuel consumption theorem) and in 
assumptions for TP-AGB stars (e.g., T c ff scale, bolometric corrections, detailed evolution, lifetimes) among 
the various works. The time evolution and difference between Lboi and Lk of the TP-AGB contribution 
in our models agree qualitatively with these other models but the comparison suggests that the peak at 
« 5 x 10 8 yr for t sc = 1 Myr may be too narrow and with too large an amplitude. Around 10 8 yr, our 
short burst predictions lie in the ranges computed for coeval populations by the above authors, from a few 
to ~ 15% for Lboi and < 10% to ~ 40% for Lk (with fairly large differences between different work). We 
point out that we implemented the TP-AGB to assess to first order the effects on the properties synthetized 
by STARS, and that our results of M82 would be essentially unchanged without it. Accurate modeling 



22 



would require a more careful and detailed treatment as presented by, e.g., Maraston (1998) and Mouhcine 
& Lancon (2002), as well as a better understanding of the properties and evolution of TP-AGB stars. 
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Fig. 1. — Comparison of the observed properties 
of selected regions in M 82 with model predictions. 
The curves are computed for m up = 100 M and 
four different burst decay timescales: 1 Myr (solid 
lines), 5 Myr (dashed lines), 20 Myr (dash-triple 
dot lines), and 1 Gyr (dotted lines). The other 
model parameters are given in Table 2. The hori- 
zontal bars indicate the measurements for the dif- 
ferent regions as follows. In the left-hand side pan- 
els: central 35 pc at the nucleus (empty bars la- 
beled "N"), Bl (light-shaded bars), and B2 (dark- 
shaded bars). In the right-hand side panels: the 
3D field of view and the entire starburst core 
(dark- and light-shaded bars, respectively). The 
width of each bar corresponds to the formal un- 
certainties. 
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Fig. 2. — Comparison of the neon and L^/L^ yc 
ratios for selected regions in M 82 with model pre- 
dictions. The curves correspond to m up = 100 M 
(solid lines), 50 M Q (dashed lines), 35 M Q (dash- 
dotted lines), 30 M Q (dash-triple dot lines), and 
25 M Q (dotted lines); two burst timescales are il- 
lustrated: 1 and 5 Myr (steepest and shallowest 
curves, respectively). The horizontal bars within 
each plot indicate measurements for different re- 
gions as in Figure 1. The vertical bars and thick 
lines on the right-hand side of the figure show 
the ranges of and average ratios for the sample of 
solar-metallicity starburst galaxies observed with 
750-SWS by Thornley et al. (2000). 
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Fig. 3. — Single burst evolutionary tracks in dia- 
grams of the neon, Lboi/^Lyc, and Lk/L^ yc ratios 
versus W-2.29- The data points with error bars in- 
dicate the measurements for the nucleus (filled cir- 
cles), Bf (filled triangles), B2 (filled squares), the 
3D field of view (open triangles) , and the starburst 
core (open circles) . The thick and thin curves cor- 
respond to models with m up = 100 and 50 M , 
respectively. Two burst decay timescales are plot- 
ted for each m up : 1 and 5 Myr (solid and dotted 
line pairs, respectively). The open squares along 
each curve indicate burst ages separated by loga- 
rithmic intervals of Alog(ib [yr]) = 0.1 dex, with 
selected ages (in Myr) labeled for reference. 
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Fig. 4. — Two-burst model evolutionary tracks in 
diagrams of the neon, Lboi/iLyc, and Lk/L^c 
ratios versus W2.29- The data points with error 
bars indicate the measurements for the nucleus 
(left panels), Bl and B2 (center panels), and the 
3D field of view and starburst core (right panels). 
The curves show the evolution of the properties 
for the combination of two bursts with best-fit 
time separation and relative intensities for each 
region (§ 4.4), each burst with t sc = 1 Myr and 
m up = 100 M Q . Symbols and lines are as follows: 
filled circle and dotted line for the nucleus, filled 
triangle and short-dashed line for Bl, filled square 
and long-dashed line for B2, open triangle and thin 
solid line for the 3D field of view, and open circle 
and thick solid line for the starburst core. The 
open squares and labels along each curve indicate 
the time in Myr since the onset of the older burst, 
in logarithmic intervals of Alog(tb [yr]) = 0.1 dcx. 




Fig. 5. — Comparison of the CO bandheads EWs 
of selected regions in M 82 with model predictions. 
The data and model curves are the same as in 
Figure 1. A larger range of burst ages is shown 
to compare the EWs during the phases when red 
supcrgiants dominate the near-IR continuum and 
when intermediate-mass stars become important 
contributors, in particular TP-AGB stars. The 
vertical bars and thick lines on the right-hand 
side of the figure indicate the ranges of and aver- 
age EWs observed for normal stellar populations 
characteristic of elliptical galaxies and bulges of 
spiral galaxies (Oliva et al. 1995). The sharp in- 
crease in synthetized EWs for t sc < 20 Myr around 
500 Myr is most likely exaggerated due to the syn- 
thesis technique employed (see appendix A). Cor- 
responding solutions are thus not considered valid 
and are also not supported by the mass-to- JT-band 
light ratio. 
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Fig. 6. — Starburst models for selected regions in 
M82. The curves in each plot represent the evo- 
lution of the properties for the two-burst model 
for each region (§ 4.4) as a function of the time 
elapsed since the onset of the oldest burst: neon 
ratio (grey solid lines labeled "Ne"), Lboi (black 
dotted lines), L^ yc (black dashed lines), Lk (black 
solid lines), and W\^2 and W2.29 (grey and black 
dash-dotted lines, respectively). The ages for the 
bursts ("old" and "young") are given in each dia- 
gram. The curves are normalized to the observed 
values; ideally, they should all meet at unity (hori- 
zontal line) at the appropriate age for the old burst 
(vertical line). An m up = 100 M Q and t sc = 1 Myr 
were adopted for each burst. 




Burst age [yr] 



Fig. 7. — Effects of different model assumptions 
and parameters on the predictions of selected 
properties. For all models shown, m up = 100 M 
and t sc = 1 Myr. The black solid curves repre- 
sent the models applied to M82, computed as- 
suming a homogeneous cluster population. The 
other curves show the effects of accounting for a 
heterogeneous cluster population following a plau- 
sible luminosity function (grey solid lines), and 
of adopting the IMF of Miller & Scalo (1979, 
dashed lines) or Eisenhauer et al. (1998, dash- 
dotted lines). Additional curves for the neon ra- 
tio were computed for a gas and dust composition 
typical of the Orion nebula (dashed-triple dot line) 
and for log U e g = — 2 dex and —2.5 dex (upper and 
lower dotted lines, respectively). The properties of 
the starburst core of M 82, representative of those 34 
of the other selected regions, are indicated by the 
horizontal bars. 




Fig. 8. — Comparison of various IMFs, repre- 
sented by different lines as labeled in the plot. 
The IMFs are normalized to unity at I M Q . The 
shaded area indicates the range of slopes deter- 
mined in young clusters and OB associations in 
the Milky Way and in the Large Magellanic Cloud 
(Hunter ct al. 1997 and references therein; Brandl 
ct al. 1996; Eisenhauer et al. 1998; Massey & 
Hunter 1998; Scalo 1998). 
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Fig. 9. — Predicted mass involved in the star- 
burst as a function of m\ ow . The curves show 
the burst masses for the nucleus and the star- 
burst core of M82 (central 35 and 500 pc) com- 
puted from the respective two-burst models nor- 
malized such that the predicted luminosities (Lk 
and I/Lyc) equal the observed luminosities for the 
derived burst ages (see Tables 3 and 4). Horizon- 
tal lines mark various fractions of the total stellar 
mass determined for each region (7.9 x 10 7 M 
and 6.1 x 10 8 M for the central 35 pc and the 
starburst core, respectively). The top and middle 
panels show, for each region, the effects of adopt- 
ing different IMFs, represented by the different 
curves as labeled in the plots. The bottom panel 
illustrates results obtained for the starburst core 
with a Salpeter (1955) IMF but assuming different 
values for the if -band luminosities (solid curves; 
"FS01": this work, "M93": McLeod et al. 1993, 
"S97": Satyapal et al. 1997; see § 4.6); the dashed 
curve show computations assuming the Rieke et 
al. (1993) IMF and McLeod et al. (1993) L K . 
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Fig. 10. — Properties of individual 1" x 1" pix- 
els within the 3D field of view. The diagrams 
are as in Figure 3, with model curves computed 
for a single burst with m up = 100 M Q , and 
t sc = 1 and 5 Myr (solid and dashed lines, respec- 
tively); open squares indicate burst ages separated 
by Alog(ib [yr]) =0.1 dex, with selected ages (in 
Myr) labeled for reference. The size of the data 
points is proportional to the intrinsic stellar Lk- 
Typical uncertainties are shown by the error bars 
in each plot. 



37 




Fig. 11. — Model results for individual pixels in 
the 3D field of view. The top panels (a and b) 
show the derived ages for the young and old bursts 
while the bottom panels (c and d) show the ini- 
tial star formation rates (Ro)- The cross marks 
the position of the nucleus and square boxes indi- 
cate selected regions modeled in § 4: the central 
35 pc at the nucleus and regions Bl and B2. The 
diagonal lines indicate the 3" -wide slit along the 
galactic plane of M 82 used to extract the profiles 
shown in Figure 12. 
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Fig. 12. — Variations of the burst age and strength 
with projected radius from the nucleus of M82. 
The radial profiles are obtained from the results 
of individual pixels in the 3" -wide slit along the 
galactic plane shown in Figure 11. The results for 
the young and for the old bursts are represented by 
open and filled circles, respectively. The positions 
of the nucleus and of the regions Bl and B2 are 
indicated. 
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Fig. 13. — Global star formation history within 
the 3D field of view in M 82. Top panel: integrated 
star formation rate derived from the spatially de- 
tailed modeling and double Gaussian profile which 
provides a good approximation thereof (solid line 
with filled dots and dashed line, respectively; see 
§ 5.2). The best-fit two-burst model for the global 
properties of the 3D field of view is plotted for 
comparison (from Table 4; dotted line). Middle 
panel: initial star formation rate surface density, 
obtained by dividing the Rq versus burst age curve 
of the top panel by the total area of the pixels con- 
tributing in each age bin. Bottom panel: model re- 
sults for the individual pixels, with open and filled 
circles representing the young and old bursts, re- 
spectively Typical error bars are indicated at the 
bottom of the plot. 
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Fig. 14. — Sketch of the succession of events re- 
lated to starburst activity in M82 depicted in 
§ 6.2. The galaxy is illustrated as viewed from 
above the plane except for the bottom panel, 
where it is inclined to show the starburst wind. 
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Fig. 15. — Selected evolutionary tracks in 
the theoretical Hertzsprung-Russell diagram for 
intermediate-mass stars. Black lines represent the 
Geneva tracks. Grey lines show our extensions 
to the TP-AGB phase computed from the pre- 
scriptions of Bedijn (1988); the star symbols cor- 
respond to the two evolutionary points associated 
with the different pulsation modes considered (ap- 
pendix A). Total TP-AGB lifetimes are given in 
parenthesis. 
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Fig. 16 — Effects of TP-AGB stars in our star- 
burst models. The left panel shows the relative 
contribution of the TP-AGB phase to the total 
if -band and bolometric luminosities (black and 
grey lines, respectively). The right panel com- 
pares the model predictions for the EW of the 
CO bandhead at 2.29 ^m accounting for TP-AGB 
stars (black lines) and excluding TP-AGB stars 
(grey lines). Computations are shown for burst 
decay timescales of t sc = 1 Myr and 1 Gyr (solid 
and dotted lines, respectively). The effects of TP- 
AGB stars at ages ~ 500 Myr for t sc = 1 Myr 
are likely exaggerated due to the conventional 
synthesis technique employed in our models (ap- 
pendix A). 
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Table 1 

Observational constraints for selected regions in M 82 



Property 


Units 


Nucleus 


Bl 


B2 


3D field 


Starburst core 


L K 


10 


8 T 

Lq 


0.56 ±0.23 


0.074 ± 0.024 


0.15 ±0.05 


2.7 ±0.8 


13 ±4 


I/Lyc 


10 


8 L 


0.67 ±0.25 


1.5 ±0.4 


2.0 ±0.9 


27 ± 13 


77 ±23 


ibol 


10 


8 T 

Lq 


18 ±9 


8.9 ±3.1 


12 ±4 


220 ± 90 


660 ± 120 


rOB 


10 


8 T 

Lq 


6.6 ±2.0 


7.4 ±2.2 


9.2 ±2.7 


170 ± 60 


390 ± 40 


M* 


10 ! 




n 7q+o- 22 








6 1+ 2 - 7 


a 


10" 


V" 1 


0.25 


0.25 


0.58 


13 


6 






A 


5.6 ±0.3 


3.4 ±0.3 


4.6 ±0.3 


4.8 ±0.3 


4.8 ±0.3 


W2.29 




A 


15.2 ± 1.2 


11.2 ± 1.4 


14.4 ± 1.4 


14.4 ± 1.4 


14.4 ± 1.4 


LK / L-Ly C 






0.84 ±0.46 


0.049 ±0.021 


0.075 ± 0.042 


0.10 ±0.06 


0.17 ±0.07 


I/bol /-I/Lyc 






27 ± 17 


5.9 ±2.6 


6.0 ±3.4 


8.1 ±5.1 


8.6 ±3.0 


^bol/^Lyc 






9.9 ±4.7 


4.9 ± 2.0 


4.6 ±2.5 


6.3 ±3.8 


5.1 ± 1.6 


M*/L K 


M c . 


)/Lq 


1.4 ±0.7 








0.47 ±0.25 


10 12 ^SN/Lbol a 


yr 




1.4 


2.8 


4.8 


5.9 


0.91 


[Ne III] /[Nell] b 






0.13 


0.16 


0.12 


0.16 


0.16 ±0.04 



a Uncertainties of a factor of two are adopted for ^sn and 10 12 ^sN/I/boi- 

Values for the nucleus, Bl, B2, and the 3D field of view are the inferred "equivalent ratios" as described 
in § 2. The uncertainties are estimated to be a factor of two. 
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Table 2 

Summary of model parameters 



Parameter Symbol Value 

IMF: dN/dm oc m" Q 



Power-law index a 2.35 (Salpeter 1955) 

Upper mass cutoff m up 25 M© — 100 M Q 

Lower mass cutoff m\ ow 0.1 M Q — 5 M© 





Star formation history 




Burst timescale 


t sc 1 Myr - 


1 Gyr 


Burst age 


tb 1 Myr - 


100 Myr 


Burst strength a 


Ro 






Stellar properties 




Metallicity 


Solar 




Mass-loss rate 


Normal 






Nebular properties 




Electron density 


n c 300 cm" 


3 



Ionization parameter* 3 \ogU c s — 2.3 dex 

Inner radius b R c g 25 pc 

Gas-phase abundances ... Solar 

Dust within nebulae ... Neglected 



a Characterized by the initial star formation rate. 

b Effective values representing the local nebular condi- 
tions in the idealized thin gas shell geometry (see § 3.2). 
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Table 3 

Starburst models for the nucleus and regions Bl and B2 



Nucleus 



Parameter 




Units 


Younj 


; burst 


Old burst 


t b 




Myr 


4.7 ±0.5 


9 6 +u - B 

M -°-0.2 




Ro 




M yr- 1 


0.52 


+0.30 
-0.25 


6.6 ± 1.5 


Property 




Units 


Observed 


Young burst 


Old burst 


Total 


Lk 




10° L 


0.56 ± 0.23 


0.011 


0.55 


0.56 


I/Lyc 




10 8 L 


0.67 ±0.25 


0.48 


0.16 


0.64 


ibol 




10 8 L 


18 ±9 


6.8 


27 


34 


M* 




10 8 M 


0.79±°; 2 2 ? 


0.0052 


0.066 


0.071 


Wl.62 




A 


5.6 ±0.3 


0.0 


5.6 


5.5 


W2.29 




A 


15.2 ± 1.2 


0.0 


15.2 


14.9 


a 




10- 2 yr _1 


0.25 


0.038 


0.45 


0.49 


[Ne III] /[No 


II] a 




0.13 


0.25 


0.0038 


0.21 


Bl 


Parameter 




Units 


Younf 


; burst 


Old burst 


t b 




Myr 


3.9 


+0.4 
-0.6 


7.8i^ 




-Ro 




M Q yr- 1 


0.84 ± 0.26 


96+ ' 27 


Property 




Units 


Observed 


Young burst 


Old burst 


Total 


Lk 




10 s L 


0.074 ± 0.024 


0.013 


0.061 


0.074 


L-Lyc 




10 8 L Q 


1.5 ±0.4 


1.4 


0.083 


1.5 


Lbo\ 




10 8 L Q 


8.9 ±3.1 


14 


5.7 


20 


M* 




10 8 M 




0.0082 


0.0096 


0.018 


Wl.62 




A 


3.4 ±0.3 


0.0 


3.5 


2.9 


W2.29 




A 


11.2 ± 1.4 


0.0 


10.9 


9.0 


a 




htV- 1 


0.25 


0.040 


0.078 


0.12 


[Ne III]/[Ne 


II] a 




0.16 


0.45 


0.017 


0.42 


B2 


Parameter 




Units 


Younf 


; burst 


Old burst 


t h 
Ro 




Myr 
M Q yr" 1 


4.0 
1.2 


+U.4 
-0.5 
+0.5 
-0.4 


8.9 ±0.2 

1 7+ ' 4 
L - 1 -0.5 


Property 




Units 


Observed 


Young burst 


Old burst 


Total 


Lk 




10 s L 


0.15 ±0.05 


0.019 


0.13 


0.15 


LLyc 




10 8 L 


2.0 ±0.9 


1.8 


0.067 


1.9 


Lbo\ 




10 8 L 


12 ±4 


19 


8.0 


27 


M* 




10 8 M 




0.012 


0.017 


0.029 


Wl.62 




A 


4.6 ±0.3 


0.0 


4.9 


4.2 


W2.29 




A 


14.4 ±1.4 


0.0 


13.8 


12.0 


Z^SN a 




htV- 1 


0.58 


0.062 


0.13 


0.19 


[Ne III]/[Ne 


II] a 




0.12 


0.41 


0.0067 


0.40 



a The adopted uncertainties for the constraints on the neon ratio and i^sn are a 
factor of two. 
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Table 4 

starburst models for the 3d field of view and the starburst core 



3D field of view 



Parameter 




Units 


Young burst 


Old burst 


tb 
-Ro 




Myr 
M yr- 1 


4.1^ 


9.0 ±0.2 

3111 


Property 




Units 


Observed 


Young burst 


Old burst 


Total 


L K 




10 s L Q 


2.7 ±0.8 


0.30 


2.4 


2.7 


Lhyc 




10 8 L Q 


27 ± 13 


25 


1.1 


26 


ibol 




10 8 L 


220 ± 90 


280 


140 


420 


M* 




10 s M 




0.18 


0.31 


0.49 


Wl.62 




A 


4.8 ±0.3 


0.0 


5.0 


4.4 


W2.29 




A 


14.4 ± 1.4 


0.0 


14.0 


12.4 


^SN a 






13 


0.99 


2.3 


3.3 


[Ne III]/[Nc 


II] a 




0.16 


0.38 


0.0061 


0.37 


Starburst core 


Parameter 




Units 


Younj 


j burst 


Old burst 


*b 




Myr 


3.9 


+U.2 
-0.3 


9.0 ±0.2 


Ro 




M yr" 1 


43 ± 11 


160 ± 30 


Property 




Units 


Observed 


Young burst 


Old burst 


Total 


L K 




10 s L Q 


13 ±4 


0.66 


12 


13 


Lhyc 




10 8 L 


77 ±23 


69 


5.8 


75 


Lboi 




10 8 L Q 


660 ± 120 


710 


740 


1500 


M* 




10 8 M 


6 1+ 2 ' 7 

°- 1 -2.5 


0.42 


1.6 


2.0 


Wl.62 




A 


4.8 ±0.3 


0.0 


5.0 


4.7 


W2.29 




A 


14.4 ± 1.4 


0.0 


14.0 


13.3 


a 




ictV- 1 


6 


2.1 


12 


14 


[Ne III]/[Ne 


II] 




0.16 ±0.04 


0.45 


0.0061 


0.42 



a The adopted uncertainties for the constraints on the neon ratio for the 3D field 
of view and on i^sn for both regions are a factor of two. 
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Table 5 

Parameters for various IMFs a 



Source 


an 


mi 
(M ) 


OC2 


m 2 
(M ) 


Q!3 


Salpeter 1955 


2.35 










Miller & Scab 1979 


1.4 


1 


2.5 


10 


3.3 


Rieke et al. 1993 


1.25 


3 


2.5 


10 


3.0 


Eisenhauer et al. 1998 


1.73 


15 


2.7 







a The IMFs are represented by a power-law diV/dm oc 
mT a . The values of ai give the power-law indices for dif- 
ferent mass ranges and the m, correspond to the inflection 
points for broken power-laws. 
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